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CHAPTER  1   Introduction 

1.1  GP  and  Traditional  HODH  Approaches 

Decision  making  la  the  process  of  selecting  a  possible 
course  of  action  from  all  the  available  alternatives.  In 
almost  all  such  problems  the  multiplicity  of  criteria  for 
judging  the  alternatives  is  pervasive.  That  is,  for  many 
such  problems,  the  decision  maker  wants  to  attain  more  than 
one  objective  or  goal  In  selecting  the  course  of  action 
while  satisfying  the  constraints  dictated  by  environment, 
processes,  and  resources.  Another  characteristic  of  these 
problem  is  that  the  objectives  are  apparently 
non-commensurable . 

Traditionally  there  are  two  approaches  for  solving  the 
multiobjectlve  problem.  One  of  them  is  to  optimize  one  of 
the  objectives  while  appending  the  other  objectives  to  a 
constraint  set  so  that  the  optimal  solution  would  satisfy 
these  objectives  at  least  up  to  a  predetermined  level.  The 
problem  Is  given  as  :  (  assume  that  all  the  objectives  are 
for  maximization  ) 

max   f^x)  (l.l) 

s.t.    g^Cx)  <  0  ,  j  =  1 m 

f^(x)  >  a   ,  Jt»  1,  ....  k,  and  i*  1 
where  a^  is  the  acceptable  predetermined  level  for  objective 


jL.      The  other  approach   13   to  minimize  a   super-objective 

function  created  by  multiplying  each  objective  function  with 

a  suitable  weight  and  then  by  adding  them   together.     This 

approach  leads  to  the  solution  of  the  following  problem  : 

IC 

maX    I     wifi(ii>  (1.2) 

1  =  1 

s.t.      gj(x)   <  0.    j    =   1 m. 

k 
The  weights  are  usually  normalized  so  the  S  w   =1. 

1  =  1 

Both  of  the  above  approaches  are  ad  hoc  at  best.    Often 

they  lead  to  a  solution  which  may  not  be  the   beat  or  most 

satisfactory.   Because  of   the   incommensurability  and   the 

conflicting  nature  of  the  multiple   criteria,   the  problem 

becomes  complex  and   it   becomes   difficult   to   choose   the 

acceptable  levels,  a^'s.  in  (1.1)   that   will   result   in  a 

nonempty   constraint   set   in   the   first   attempt   for   the 

solution.     In  the   first  approach,   the   Implied   value 

trade-off  between  f  .  and  f   is  as  follow  : 


L     »,       Ijl    <    a 


Value   Trade-off    = 

U   <  a/ 

This  may  not  be  the  actual  value  structure  of  the  decision 
maker,  and  this  value  structure  is  sensitive  to  the  level   of 
a^.   For  the   second   approach,   the   major   problem   is   in 
determining  the   proper   weight,   w±,   especially   when   the 
objectives  are  Incommensurable  and  incompatible.     Multiple 


Objective  Decision  Making  (  MODM  )  methods  are  the  result  of 

the  desire  to  eliminate  the  above  difficulties. 

Goal   Programming   (GP)   was   originally   formulated   by 

Charnes  and  Cooper  [CI]  in  1952,  which  tried  to   Incorporate 

the  above  two  approaches  without  their  difficulties.   The  GP 

method  requires  the  decision  maker  to  set  the  goals  for  each 

objective  that  he/she  wishes  to  attain,   and   introduce   the 

positive  and  negative  deviational   (   surplus   and   slack   ) 

variables  to  avoid  the  situation  of  empty  constraint  set.   A 

preferred   solution  is   then  defined   as   the  one   which 

minimizes  the  deviations   from  the   set  of   goals.     The 

formulation  is  given  as  s 

k 

min  [  I     widI 
1  =  1 

s.t.   g.(x)  <  o,  j  =  i,  2.    ....   m 

fi(-)  +  di  "  di  =  ai'  L   =  l»  2 k 

To  avoid  the  difficulty  of  specifying  the  weights,  Ijiri 

[13],  in  1965,  introduced  the  definition  of  preemptive 
priority  to  GP.  The  decision  maker  does  not  need  to  give 
numerical  weights  for  the  objectives;  he/3he  only  needs  to 
give  an  ordinal  ranking  of  them.   This   version   of   GP   was 

very  popular  in  solving  the  MODM  problem.   It  seemed  that  it 
solved   the   difficulties   of   the   traditional   approaches. 
Actually,  it  did  not.   GP  still  has  the  same  problems  as  the 
traditional  approaches.   We  are  going  to  discuss  it  in   this 


paper. 


1.2  Research  Objective 

GP  has  been  a  popular  MODM  method  and  has  been  applied 
to  almost  every  field  in  solving  the  MODM  problems.  But 
some  inherent  shortcomings  of  GP  are  not  apparent  and 
well-known  to  its  users.  People  who  use  a  method  and 
without  knowing  its  limitations  are  dangerous.  Therefore, 
the  major  purpose  of  this  study  is  to  clarify  the 
shortcomings  of  GP.  Although  this  study  is  not  the  first 
attempt  to  point  out  the  shortcomings  of  GP,  this  study 
makes  it  more  complete  and  detailed. 

In  addition,  some  remedial  methods  are  suggested  to 
enhance  the  adaptability  of  GP.  Also,  efforts  are  made  to 
propose  a  general  decision  rule  to  help  the  user  to  identify 
whether  GP  is  the  appropriate  method.  If  GP  is  not  the 
appropriate  method,  then  the  other  MODM  methods  are 
recommended . 


1.3  Thesis  Structure 

This  report  discusses  the  formulation,  algorithms, 
shortcomings,  and  remedial  methods  of  GP.  Chapter  1  briefly 
introduces  the  relationship  of  GP  and  traditional  approaches 
to  MODM  problems  and  the  purpose  of  the  study.  Chapter  2 
reviews  the  linear  programming,  the  MODM  problem 
formulation,  classification  of  MODM  methods,  and  the 
different  GP  formulations  and  algorithms.  Chapter  3 
describes  the  shortcomings  of  GP  and  illustrates  these 
shortcomings  by  numerical  examples.  Chapter  4  presents  the 
suggested  remedial  methods  for  GP  and  the  decision  rule  for 
recommending  the  other  MODM  methods  in  different  conditions. 
Finally  the  conclusions  are  given  in  Chapter  5. 


CHAPTER  2  Goal  Programing  and  HODM 

Goal  Programming  is  an  extension  of  linear  programming 
originally  formulated  by  A.  Charnes  and  W.  W.  Cooper  in 
1952  [CI].  It  is  also  one  of  the  MODM  (  Multiple  Objective 
Decision  Making  )  techniques.  Therefore,  the  purpose  of 
this  chapter  is  mainly  to  review  the  general  formulation  of 
LP  and  MODM  problems  and  different  types,  formulations.  and 
algorithms  of  GP. 

2.1  Review  of  Linear  Programming 

Linear  Programming  is  one  of  the  most  popular  and  well 
known  mathematical  programming  methods  which  are  U3ed  to 
solve  optimization  problems.  The  general  form  of  the  LP 
model  with  m  constraints  and  n  variables  (  n  >  ra  )  is  : 

Max/Min   Z  =  ex   +  c,x,  +  ...  +  c  x 
11    2  2  n  n 


s.t. 


*11*1  +  al2X2  +  ••••  +  alnxn  *  bl  <2-D 

'21*1  +  a22X2  +  ••••  +  a2nXn  =  b2 


a_,x   +  a  „x_  +  ....  +  a   x   =  b 

">1  1     m2  2  mn  n    Dm 

x,  b  >  0 


Note   that   although   the    constraints    are    written   as 
equalities,  mathematical  programming  does   admit   inequality 


constraints  as  well.   The  main  properties  of   the   LP   model 
are  : 

(1)  Only  single  objective 

(2)  All  variables  are  nonnegative 

(3)  Both  the  objective  function  and  constraint  functions 
are  linear 

(4)  There  always  exists  one  and  only  one      "best" 
solution,  except  in  special  Instances  (  such  as  tie, 
degeneracy,  unbounded,  or  alternate  optima  ). 

The  simplex  method,  which  was  developed  by  George  B. 
Dantzig  in  1947.  is  the  general  procedure  for  solving  LP 
problems.  It  has  been  widely  and  successfully  used  for 
solving  real-world  problems.  In  addition,  a  lot  of 
convenient  computer  codes  were  developed  and  available. 
Coverage  of  LP  solution  procedures  can  be  found  in  most  of 
the  operations  research  books  such  as  Hlller  and  Liberman 
[H5],  and  Wagner  [Wl]. 

LP  is  the  basic  tool  in  solving  the  linear 
multlobjectlve  programming  problems.  For  the  purpose  of 
describing  the  relationship  of  LP  with  MODM  and  GP  problems, 
the  same  numerical  example  (Hardee's  production  scheduling 
problem  )  is  used  in  this  chapter. 

Numerical  Example  2.1  (  Hardee's  Scheduling  Problem  ) 

The  Hardee  Toy  Company  makes  two  kinds  of  toy  dolls. 
Doll  A  is  a  high  quality  toy  and  Doll  B  Is  of  lower  quality. 


The  respective  profits  are  90.40  and  90.30  per  doll.  Each 
doll  of  type  A  requires  twice  as  much  time  as  a  doll  of  type 
B,  and  if  all  dolls  were  of  type  B,  the  company  could  make 
500  per  day.  The  supply  of  material  is  sufficient  for  only 
400  dolls  per  day  (  both  A  and  B  combined  ).  Assuming  that 
all  type  A  and  type  B  dolls  the  factory  makes  can  be  sold, 
how  should  the  manager  schedule  production  to  achieve  the 
most  profit  ? 

Thi3  problem  is  a  typical  LP  problem  which  is  formulated 
as  below  : 

max    f^x)  =  0.4x   +  0.3x 


subject  to 


xx  +  x2  <  400 


2x^    +  x   <  500 


X,,    x2  >  0 


where  x^^  and  x2  are  the  numbers  produced  of  Doll  A  and   Doll 
B  respectively.   The  optimum  solution  is. 


x1  =  100      x   ■  300 


and  the  maximum  profit  =  f  (x)  =  3130 


2.2  Hultiobjective  Programing 

2.2.1  Formulation  of  HODH  Problems 

Multiobjective  programming  deals  with  optimization 
problems  with  two  or  more  objective  functions  subject  to  a 
constraint  set.  The  multiobjective  programming  problem 
differs  from  the  LP  only  in  the  expression  of  their 
respective  objective  functions.  The  general  multiple 
objective  decision  making  problem  will  be  defined 
mathematically  as  follows  i 


max/min 
subject  to 


ff^x).  f2(x) fk(x)  1  (2.2) 


x  e  X 


where  X  =  {  x  |  g^x)  {>.=.<}  0.  i  =  1,2,..., a  j 

f  (x)  :  the  benefit  objective  for  maximization,  j  e  J 
f  (x)  :  the  cost  objective  for  minimization,  i  m   I 


Multiobjective  problems  arise  in  the  design,  modeling, 
and  planning  of  many  complex  resource  allocation  systems. 
For  example, in  the  area  of  water  resource  allocation,  the 
operation  of  a  multipurpose  reservoir  may  call  for 
delivering  irrigation  water  and  supplying  electric  power  to 
a  nearby  community,  while  still  trying  to  maintain  minimum 
water  levels  in   the   reservoir   itself   and   downstream   to 


accommodate  environmental  and  recreational  goals.  These 
goals  may  conflict  in  nature  and  be  incommensurable  and  then 
there  is  generally  no  single  optimal  solution.  Therefore, 
in  MODM  problem,  instead  of  seeking  a  single  optimal 
solution,  a  set  of  "nondominated"  solutions  is  sought. 

Nondominated  Solutions  : 

The  set  of  nondominated  solutions  is  a  subset  of  the 
feasible  region.  A  feasible  solution  in  MODM  is 
nondominated  if  there  exists  no  other  feasible  solution  that 
will  yield  an  improvement  in  one  objective  without  causing  a 
degradation  in  at  least  one  other  objective.  The  concept  of 
a  nondominated  solution  also  appears  in  the  literature  under 
the  names  of  Pareto-optimal  solution,  nonlnferlor  and 
efficient  solution. 

Preferred  Solution  : 

For  most  of  the  MODM  problems,  the  number  of  solutions  in 
the  set  of  nondominated  solutions  is  quite  large.  So,  the 
decision  maker  must  make  a  final  selection  of  the  most 
satisfactory  solution  by  using  some  other  criteria.  This 
final  solution  is  known  as  the  preferred  solution. 

The  Positive- Ideal  Solution  (PIS)  : 

In  a  MODM  problem,  the  positive-ideal  solution  is  the 
one  that  optimizes  each  objective   function   simultaneously. 
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e.g. ,  max  f  .  (x) 


xeX 


3    - 


J  and  min  £  (x) , 


I.     An   ideal 


solution  can  be  defined  as  A   =  (  f  .  f  , 


13  a  feasible  and  optimal   value   foe   the 
function.   The  PIS  payoff  table  is  given  as  : 


f   )  where  f » 


j2    objective 


fl 

f2     •   - 

•    fk 

fl 

* 

f12   •  • 

•    fik 

£2 

!21 

f2    '  ' 

•      'n 

fk 

fKl 

fk2   •  " 

* 

(diagonal  elements 
are  the  PIS) 


The  Negative- Ideal  Solution  (NIS)  : 

Similarly,  the  concept  of  a  negative-ideal  solution,   or 

the  worst  possible  solution  for  MODM  is  presented.   The   NIS 

is  the  one   that    min  f  .(x),  j  e  J   and    max  f.(x),  1  «  I, 
xeX   -1  xsX   A 

simultaneously.   A  NIS  can  be  defined  as  A~  =  (  f~,    £~ ,    .... 

fk  '  where  f^  is  a  feasible  and  the  worst  value  for  the  ^ tn 

objective  function.   The  NIS  payoff  table  is  given  as  : 


fi 

f2    •  • 

•    fk 

fl 

£l 

Z12   "  • 

•    2lk 

fa 

221 

f2    '  ' 

•    22k 

fK 

Zkl 

Zk2   •  • 

•    fk 

(diagonal  elements 
are  the  NIS) 
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Numerical  Example  2.2  :  (  Continued  from  Example  2.1  ) 

Suppose  that  the  best  customer  of  the  company  wishes   to 
have  as  many  as  possible   of   type   A   doll.     The   manager 

realizes  that  two  objectives   (1)   the   maximization   of 

profit,  and  (2)  the  maximum  production  of  Doll  A  —   should 
be  considered  in  scheduling  the  production. 

This  MODM  problem  can  be  formulated  as  follows  : 
max  f,(x)  =  0.4x   +  0.3x, 

max  f2(x)  =  x1 


subject  to 


x   +  x   <  400 


2x1  +  x£  <  500 


V  X2  *  ° 


Step  0-1.  Obtain  the  PIS  and  construct  the  PIS  payoff  table. 

1)  max  f  (X)  =  0.4x   +  0.3x, 
xeX  1        Z 

This  is  a  LP  problem.   The  solution  is  : 

*    *  * 

(x^  x2)  =  (100,  300),  f  (x)  =  130 

2)  max  f,(x)  =  x 
xsX   Z        X 

This  is  again  a  LP  problem.   The  solution  is  : 
(x1,  x2)  =  (250,  0),  f*(x)  =  250 
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I   "2 
AC  0  ,400) 

(100,300) 

C(250,  0  ) 


Figure  2.1  The  decision  variable  space 
representation . 


f  (PIS) 


1    2 
AC120,  0  ) 


0. 4x, +0. 3x. 


BC130,100) 

C(100,250) 


Figure  2.2  The  objective  function  space 
representation . 
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These  results  are  summarized  in  the  PIS  payoff  table 


fl 

f2 

Xl 

X2 

fl 

* 
130 

100 

100 

300 

f2 

100 

250* 

250 

0 

Step  0-2.  Obtain  the  NIS  and  construct  the  NIS  payoff  table. 


1)  min  f  (x)  =  0.4x   +  0.3x, 
xeX  1        L 


The  solution  is 


(x1#  x2)  =  (  0  .    0  ),  f^x)  =  0 


2)  min  f,(x)  =  x. 
xeX  L  1 

The  solution  is  : 

(xl'  x2)  =  (  °  '  °  )#  f2(-)  =  ° 
These  results  are  summarized  in  the  NIS  payoff  table 


fl 

f2 

xl 

X2 

fl 

0~ 

0 

0 

0 

f2 

0~120 

o" 

0 

0~400 

The  results  are  illustrated  in  Figure  2.1  and  2.2.  The 
set  of  nondominated  solution  is  on  BC  and  the  preferred 
solution  should  be  one  point  on  BC. 
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2.2.2  Classification  of  MODM  Methods 

To  solve  the  MODM  problem,  a  lot  of  methods  have  been 
created.  Each  method  has  different  types  of  assumptions, 
and  different  information  requirements.  In  the  process  of 
decision  making,  some  preference  Information  articulated  by 
the  decision  maker  may  be  required  in  order  to  choose  the 
preferred  solution  out  of  a  set  of  nondominated  solutions. 
Hwang  and  Masud  [H6]  classified  solution  techniques  for  MODM 
problems  according  to  the  timing  of  the  requirement  for 
preference  information  versus  optimization.  Figure  2.3 
depicts  that  scheme  of  classification  in  detail, 
classification  in  detail. 

The  first  class  requires  no  information  from  the  DM  once 
the  problem  constraints  and  objectives  have  been  defined. 
The  analyst  makes  assumptions  about  the  DM's  preference  and 
presents  solution(s)  to  the  DM. 

The  second  class  assumes  the  DM  has,  consciously  or 
subconsciously,  a  set  of  goals  to  achieve  and  the  goals  set 
are  given  to  the  analyst  before  the  formulation  of  a 
mathematical  model.   GP  belongs  to  this   class. 

The  third  class,  also  known  as  interactive  methods, 
requires  greater  DM  involvement  in  the  solution  process. 
The  interaction  takes  place  through  a  DM-analyst  or 
DM-computer  dialogue  at  each  iteration.  The  tradeoff  or 
preference  information  given  by  the  DM  after  each  iteration 
is  used  for  determining  a  new  solution.  The  DM,  hence, 
actually  gains  insights  into  the  problem. 
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The  last  class  of  methods  determines  a  subset  of  the 
complete  set  of  nondominated  solutions  to  the  MODM  problem. 
In  doing  so,  it  deals  strictly  with  the  physical  constraints 
and  makes  no  attempt  to  consider  the  preference  of  the  DM. 
The  desired  outcome  is  to  narrow  down  the  possible  solutions 
and  make   the  selection  of  the  preferred  solution  easier. 

It  is  worth  noting  that  none  of  the  methods  are  superior 
to  the  others,  because  one  method's  advantage  is  another's 
shortcoming  and  vice  versa.  Therefore,  to  choose  a  MODM 
method  is  a   MODM  problem  itself. 
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2.3  Goal  Programing  (GP) 

2.3.1  Literature  Survey  of  Goal  Fi  iiji  ii—l  iiij 

The  birth  of  GP  was  in  1952  when  Charnes  et  al.  [CI] 
obtained  "constrained  regression  estimates"  for  a  problem 
involving  the  development  of  an  executive  compensation 
formula.  The  term  "goal  programming"  was  first  coined  by 
Charnes  and  Cooper  [C2]  in  their  1961  classic  text.  In 
essence,  they  proposed  a  model  and  approach  for  dealing  with 
certain  LP  problems  In  which  the  conflicting  goals  of 
management  were  inappropriately  forced  to  be  included  as  the 
constraints  of  LP  model.  Since  then,  researchers  began  to 
use  GP  for  raultiobjective  or  raultigoal  problems.  In  the 
last  15  years,  there  has  been  a  proliferation  of  articles 
pertaining  to  the  application  and  theory  of  GP.  A  major 
force  in  this  popularity  was  Lee,  who  authored  the  first 
book  devoted  exclusively  to  GP  [LI].  An  earlier  book  by 
Ijlrl  [13]  applied  GP  to  accounting  problems.  A  subsequent 
text  by  Ignizo  [II]  contributed  further  to  the  popularity  of 
GP. 

Previous  surveys  of  GP  have  been  conducted  by  Kornbluth 
[K2],  Charnes  and  Cooper  [C3],  Ignizlo  [12],  Lin  [L2], 
Hannan  [H3],  and  Zanakis  and  Gupta  [Zl].  Kornbluth's  survey 
[K2],  written  in  1973,  is  somewhat  outdated  because  there 
have  been  many  theoretical  advances  in  GP  in  the  intervening 
years.  The  Charnes  and  Cooper  survey  [C3],  conducted  in 
1977,  contains  primarily  a  survey   of   their   work   on   such 
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extensions  as  explicit  methods  in  GP,  goal  interval 
programming,  and  fractional  GP.  Ignizio's  study  [12]  is 
primarily  to  readers  not  familiar  with  GP .  It  presents  a 
historical  sketch,  lists  applications,  describes  codes,  and 
discusses  areas  for  futher  research.  Hannan's  survey  [H3] 
summarized  methodological  contributions  to  the  GP  of  the 
period  of  1973-1982  and  made  comment  about  the  merits  and 
shortcomings  of  GP  relative  to  other  multiple  objective 
decision  making  techniques.  Lin  [L2]  has  provided  a 
bibliography  containing  84  applications  of  GP  to  accounting, 
finance,  operations  management,  marketing,  and  manpower 
planning.  Zanakls  and  Gupta  [Zl],  in  1985,  further  compiled 
and  classified  240  GP  articles  according  to  techniques  and 
application  areas. 

Despite  this  rosy  picture,  GP  has  been  a  controversial 
MODM  method.  The  formulations  of  GP  hide  some  major 
shortcomings  in  solving  the  MODM  problems.  These 
shortcomings  have  been  discussed  by  some  researchers  and 
users  of  GP,  such  as  Hannan  [H2,H3],  Zeleny  [Z2,Z3],  Cohon 
and  Marks  [C5],  Harrald  [H4],  Sherali  [S2],  Sherali  and 
Soyster  [S3],  Gass  [Gl],  and  Alvord  [Al].  The  most 
recurrent  criticisms  relate  to  the  use  of  preemptive 
priorities  and  the  priori  definition  of  goals. 


2.3.Z  Formulation  of  Goal  Programming  Problems 

The  general  mathematical  representation  of  the  multiple 
objective  decision  making  problem  is  as  stated  in  section 
2.2  : 

max/min   [f^x),  f2(x),  fk(-}  1  (2-2) 

subject  to 

x  m  X 

where  X  =  |  x  |  g^x)  {>,  =  ,<}  0,  1  =  1,2 m  \ 

In  the  GP  approach  of  solving  the  problem  as  posed  by 
(2.2),  the  decision  maker  is  required  to  indicate  his 
targets  or  goals  b  .  (  aspiration  level  )  for  each  of   the   k 

objectives  f  .  (x).    In  most   cases,   it   is   impossible   to 

satisfy  exactly  all  such  goals.   The   GP,   then,   employs   a 

minimum-distance  notion  of  best.   That  is  to  minimize  the  d 

P 

distance  from  the  goals  : 
k 


min  dp  =  1   Z  '  bj  "  f j(-)  !   j     •   P  *  1         (2.3) 


subject  to 


gi(x)  <  o,   i  =  l m 

x  >  0 


where  the  b   is  the   initial   goal   for   objective   function 


e^x). 


The  GP  usually  assumes  that  the  decision   maker   has   no 
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information  on  the  ideal  solution.  One,  then,  cannot 
guarantee  that  b.  >  f  .(x)  for  all  feasible   x.     Therefore, 

the  absolute  value  is  used  in  the  distance  metric.  To  drop 
the  absolute  value  in  (2.3),  we  need  to  Introduce  two  other 
nonnegative  variables,  called  deviational  variables,  to 
substitute  the  objective  function  in  (2.3).  Then  the  new 
formulation,  i.e.,  GP  fomulation  becomes  : 


min  [  I   (  d"  ♦  d+  )P  ] 


subject  to 


gx(x)  <  0,  i  =  l, 

f  .(x)  ♦  d"  -  a+   =  „       j  =  i. 


d~  •  d+  =  0  ,       Vj 


x,  d".  d+  >  0 


where  d  and  d  are  positive  and  negative  deviations, 
respectively,  of  objective  j  from  its  goal.  If  f  .(x)  <  b., 
that  is  under-achlevement  of  goal  b.,  then  we  have  less  than 

the  desired  amount  of  objective   j;   this   is   a     negative 

deviation,  so  d"  =  b.  -  f  .(x)  and  dt=  0.    If   f  .(x)   >   b  , 
3     3     3  ~        3  3  ~        j 

then   we   have   more   than   a   desirable   amount,   that    is 

over -achievement   of   goal  j,  so  that   d+  =  -fb   -  f  fx)l 

3       j     j  -  ' 

f,(x)  -  b  and  d   =  0.   The  two  deviations  dt   and   d"   will 
3         3       3  3  3 
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never  nonzero  for  the  same  goal  since  total  deviations 
are  being  minimized  in  (2.3).  The  metric  in  (2.3)  with  p  = 
1  13  usually  used,  although  other  metrics  may  be  employed. 

Since  the  GP  is  based  on  the  "satisf icing"  philosophy 
espoused  by  Simon  [S4],  it  accepts  the  "good  enough" 
solution  rather  than  continues  to  search  for  the  optimal 
solution  (  preferred  solution  ).  That  is,  it  takes  "at 
least  b   (  i.e.  b  or  more  )"  achievement   for   the  benefit 

objectives,  "at  most  b   (  i.e.  b.  or  less  )"  achievement   for 

cost  objectives,  and   "exactly  b."  for  the  special   type   of 

objectives  such  as  balancing  the  budget.  These  achievements 
can  be  made  by  employing  over  and  under  achievement 
variables  in  the  objective  function  of  the  GP  formulation. 
Table  2.1  summarizes  the  approach  taken  to  accomplish  this 
desire  : 


Table  2.1   Model  Formulation  of  GP 


Goal  type 


*3(4)  *  b. 


Goal  Constraint 


f  .(x)  +  d  .  -  d  .  =  b. 

3  -      3     3     3 


Oeviational 
variables  to 
be  minimized 


d  + 
j 


*j<*>  *  b. 


fj(i)  ♦  «,  -  d.  =  b. 


fjCx)  -  b. 


f  .(x)  +  d  .  -  d  .  =  b. 
3  -      3     3     i 


d~+  d  + 
i         3 


The   decision  makers   may   also   wish   to   weigh     the 
deviations   to   reflect   the   relative   importance   of    the 
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objectives  and  the  relative  significance  of  positive  and 
negative  deviations.  This  can  be  accomplished  by  ordinal  or 
cardinal  weights.  The  GP  with  ordinal  (  more  precisely, 
preemptive  )  weights  13  called  the  preemptive  GP;  the  GP 
with  cardinal  weight  is  called  the  nonpreemptive  (  or 
archimedean  )  GP.  The  algorithm  of  these  two  types  of  GP 
will  be  discussed  in  the  next  two  sections. 

The  full  expression  of  the  general  linear  GP  model  may  be 
presented  as  : 

k 
Min  [  J     wr  ar(d)  ] 
r=l 

s.t. 

x  e  X 

fjCx)  +  d~  +  <j+  >  b^.  j  m  J 
fiCx)  +  d~  +  d*  =  b1.  1  e  I 
fp(x)  +  d"  ♦  d*  --  bp.       p  e  P 

where   f..(x)  •■  the  benefit  criteria,  j  e  J 

fi(x)  :  the  cost  criteria,  lei 

fp<£)  '  th«  two-sided  goal,  p  a  p 
d)  : 

ar(£)  =  dr  '   if  r  s  J 
ar(d)  =  d*  ,   if  r  m   I 


a„(d)  =  d   +  d+  ,   if  r  <=  P 

r  —      r     r 

the  relative  importance  of  goals  (  can   be 
either  preemptive  or  nonpreemptive  weight  ) 
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2.3.3  Preemptive  Goal  Programming 

In  1965,  Ijiri  [13]  presented  a  definition  of  preemptive 
priority  levels  so  as  to  treat  goals  according  to  their 
perceived  Importance.  The  goal  satisf icing  sequence  of 
preemptive  GP  is  based  on  the  preemptive  order.  That  Is, 
goal(s)  which  is  assigned  to  the  fir3t  priority  is  satisfied 
first,  and  then  the  second  priority  goal(s)  is  considered 
only  after  the  first  priority  goal(s)  is  fully  advanced.  In 
other  words,  the  solution  procedure  finds  a  subset  of  the 
feasible  region  which  satisfies  the  first  priority  goal(s), 
then  within  this  set  it  finds  the  reduced  set  which  is 
satisfactory  to  the  second  priority  goal(s)  and  so  forth. 
Accordingly,  the  goal  achievement  process  moves  to  the  next 
priority  when  it  begins  to  degrade  the  higher  priority  goal 
achievements.  The  solution  process  continues  until  the  last 
priority  goal(s),  is  considered.  The  preemptive  GP 
formulation  with     priority  levels  is  given  by 

min  a  =  [p^Cd",  d+),  P2a2(d",  d+) P^Cif,  d+)  1 


s.t.     x  e  X 


(2.4) 


f^x)  +  d~  -  d+  =  bj    ,      j  =  1 k 

dT  .  dt  =  0 

1      3 

dT  ,  d+  >  0 

3  3 

where  a  is  the  achievement  function  for   all   priorltl 
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^(d  ,  d  )  is  an  achievement  function  for  the   i     priority 


level,  i  =  1,2 


■  .    JL.      The   P,'3   are   preemptive   weights; 


that  la.  Pi  >>>  P1+1-     This   implies   that   no   number   w, 

however  large,  can  make  W*P1+1  >  P.-   Only  the  commensurable 

goals  may  be  assigned  to  the  same  priority  level.   Note  that 
the  achievement  function  is  a  vector  measurement. 

The  solution  algorithm   for  (2.4)  is  that  a  (  d~,  d+  )  is 

minimized   flrstj  let  min  a   =  a  .     Next  a  (d~.   d+)   is 

minimized,  but  in  no  circumstances  can  a   be  greater  than  a  . 

Thus  a  lower  ranking  achievement  function  cannot  be  satisfied 
to  the  detriment  of  a  higher  ranking  achievement  function. 
This  process  continues  until  a»(d~,  d+)  is  minimized.     Some 

special  computer  codes  for  linear  models  are  available.  For 
example,  in  Ignizio  [II]  and  Lee  [LI],  they  are  basically  a 
modified  simplex  algorithm  for  LP  problems.  However,  for  a 
moderate  size  problem  the  modified  simplex  algorithm  approach 
is  time  consuming,  and  needs  a  large  capacity  computer.  The 
aame  problem  can  be  solved  iteratlvely  (  see  Hwang  and  Masud 
[H6],  Dauer  and  Krueger  [Dl]  )  by  the  basic  simplex  algorithm 
for  which  many  efficient  computer  codes  are  available.  BY  an 
iterative  approach,  the  GP  model  can  be  decomposed  into  S. 
number  of  single  objective  problems.  Hwang  and  Masud  [H6] 
formulated  the  procedures  as  follows  : 
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Problem   1    :    To   find   x   =   (   x%,   x2 «3    )    so  aa   to 

rain        a,(d~.    d    ) 
3.t.      x  m  x 

«j<a)  +  a"  -  d;  =  b. 

d~.    d+    >   0 

(  for  some  j  which  are  associated  with  a   function  ) 

* 
Let  a1  be  the  optimum  achievement  for  problem   1.   i.e., 

*  -   + 

ax    -   min  a  (d  ,d  ) . 

Problem  2  :  To  find  x  so  aa  to 
min    a2(d".  d+) 
a.t.   x  e  X 

f.Cx)  +  d"  -  d+  =  b. 
J         3      J  3 

ai(l"'d+)  $   a* 
d~,  d+  >  0 
Cfor  some  j  which  are  associated  with  a1  and  a   functions) 

* 
Let  a2  be  the  solution  of  this  problem. 

Problem  3  :  To  find  x  so  aa  to 
min   a3(d~.  d+) 
a.t.   x  e  X 
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f  (x)  +  d~  -  d+  =  b. 

i  ]  3  J 


a2(f-d+)  <  a* 
d~,  d+  >  0 


(  for  some  j  which  are  associated  with  a 


a„,   and   a. 


1'  "Z'      ™~   "3 
functions  ) 


Let  a3  be  the  solution  of  this  problem. 


The  process  continues  until  we  solve  the   last   priority 
goals  which  takes  the  following  problem. 


Problem  j£  :  To  find  x  so  as  to 
min   a  (d~,  d  ) 


fjCi)  +  d~  -  d+  =  b..  .    j  =  i k 

al(d".d+)  <  a*  .  t  =  i £.x 
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Numerical  Example  2.3  s  (  Continued  from  Example  2.1  ) 

Suppose  the  beat  customer  of  the  toy  company  ordered  300 
dolls  of  type  A.  For  various  reasons  the  manager  could  not 
decline  this  order.  How  should  the  manager  switch  the 
production  schedule  to  meet  the  sudden  demand  ?  Assume  that 
the  order  conditions  are  the  same  as  Example  2.1. 

If  we  attempt  to  utilize  LP  to  maximize  profit,  the 
problem  is  formulated  as  below  : 

max  fjCx)  =  0.4x   +  0.3x 


s.t. 


xx    +    x2  <  400 


2x1  +  x2  <  500 


>  300 


V  *z  *  ° 


The  graphical  presentation  of  the  model  is  shown  in 
Figure  2.4.  Obviously, there  is  no  area  of  feasible  solution 
and  consequently  the  above  problem  is  unsolvable  by  LP.  In 
simple  words,  the  company  does  not  have  enough  capacity  to 
satisfy  the  customer's  demand. 

The  difficulty  can  be  handled  —  if  the  manager  wishes 
to  satisfy  his  desires  as  closely  as  possible  rather  than 
absolutely  satisfy  all  goals.  Let  us  assume  that  the 
manager's  priorities  are  : 

PI  :  Meet  the  best  customer's  order  as  closely  as  possible 
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Figure  2.4  Infeasible  solution  for  linear 
programming  model. 


G,:SO0  Doll  A 


Figure  2.5  Final  solnti-on  for  Example  2,3 
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(  i.e.,  produce  at  leaat  300  dolls  of  type  A  ) 
P2  s  Satisfy,  as  much  as  possible,  the  9240  profit. 

The  linear  GP  model  becomes  : 


min   a  =  (  d   ,  d,  ) 
s.t. 

f1  !     *t  +  dl  "  dl  =  30° 

f2  :  O^Xj^  +  0.3x2  +  d~  -  d*  =  240 


x   +  x   <  400 


2x1  +  x£  <  500 


x,  d",  d+  >  0 


The  above  model  is  shown  in  Figure  2.5.  The 
cross-hatched  area  is  the  feasible  region.  The  first 
priority  goal  calls  for  production  of  300  type  A  dolls  by 
minimizing  d~.   This  will  be  point  C  (  with  d"  =  50  ) .    For 

the  second  priority  goal,  we  may  proceed  from  point  C  to 
point  B,  but  this  movement  is  against  the  first  priority 
goal  (  increase  d1  ).  Hence,  point  C  has  to  be  the  final 
solution.   That  is 

x:  =  250,   x*  =  0, 

Goal  1  is  not  achieved.  (  f  (x)  =  250.  d"  =  50,  d+  =  0  ) 

Goal  2  is  not  achieved.  (  f2(x>  =  100.  d"  =  140,  d+  =  0  ) 

31 


Thi3  can  also  be  solved  by  the  iterative  algorithm   described 
above. 


Problem  1  :  (  for  the  1st  priority  level  ) 
min   a   =  d~ 

s.t.   x  s  X 

Xl  +  dl  ~  dl  =  30° 
x,  d~,  d+  >  0 

The  solution  ia  x1  =  250,  x£  =  0,  and  d"   =   50.     With 

this  result,  we  can  go  the  the  next  priority  level. 

Problem  2  :  (  for  the  2nd  priority  level  ) 
min    a2  =  d~ 


s.t. 


Xl  +  dl  ~  dl  =  300 


dx  <  50 

0.4x1  +  0.3x2  +  d~  -  d*  =  240 


The  solution  is  x1    =  250,  x   =  0,  and  d~=  140. 

The  results  of  the  above  two  sub-problem  are  the  same  as 
the  graphical  solution,  that  is  point  C  in  Figure  2.5  is  the 


final  solution. 
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2.3.4  Non-preemptive  Goal  Programing 

The  non-preemptive  (  archimedean  )  GP  utilizes  the 
cardinal  weights  to  represent  how  important  it  is  to  achieve 
each  goal  relative  to  the  importance  of  achieving  the  other 
goals.  GP  with  non-preemptive  weights  is  merely  a  LP  model 
and  can  be  solved  easily.  The  extensions  such  as  duality, 
and  sensitivity  analysis  can  be  handled  in  the  same  manner 
as  in  LP  problems.   The  formulation  is  : 

k 

min   7  w.(  d~  +  dt  ) 
j  =  l 

s.t.    x  m   X 

t^XJ    +  d"  -  d*  =  Dj  .      i    =    i.    2 k 

dj  •  d+  =  0  ,     Vj 

d~,  d+  >  0 


A  pretreatment  before  executing  the  above  program  is 
what  is  called  "goal  norming".  As  Wldehelm  [W3],  de  Kluyver 
[Kl],  and  Hannan  CH2]  have  noted,  care  should  be  taken  in 
developing  weights  associated  with  the  various  goals. 
Because  v.  represents  the  relative  importance  among  goals, 
each  deviatlonal  variable  (  i.e.  a"  ,  d+  )  should  have  a 
comparable  scale  through  the  normalization.  Normalization 
presents  a  convenient  way  to  correctly  set  the  desired 
weights.   Two  common  ways  of  scale  normalization  are   vector 
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normalization  and  linear  scale  transformation. 

Vector  normalization  can  be  achieved  by  associating   the 
same  Euclidean  norm  J|c3  [}   with   each   deviational   variable. 

The  Euclidean  norm  |c3|  =  (cjcj)1/2.  where  cJ   denotes   the 

.th 

3        goal  constraint  coefficient  vector  (  c   ,     .  c    ) 

jl        jn 

Thus  the  goal  constraint 


xx  +  2x2  +  3x3  +  d~  -  d*  =  4 
would  be  converted  to 

y±r  «a  ♦  — hr  x2  +  — ?-  x3  ♦  d"  -  «J  .    _^_ 

TH  /  14  /l4  /14 

xi  +  2x2 +  3X3 +  /IT  di  -  /IT  <  ■  4 

before  determining  the  relative  weights  to  d"  and  d+  in   the 
objective  function. 

The  linear   scale   transformation   can   be   achieved   by 
making  the  variation  between  the  positive  ideal   solution   ( 

* 
f   )  and  negative  ideal  solution  (  f"  )  equal  to  1   for   all 

objective  function  f  (x)  ,  j  =   l k.        This 

normalization  scheme  is  used  in  Hwang  and   Masud's   ISGP   II 

[H7].   The  goal  constraint  associated  with  the  jth  objective 
under  this  scheme,  is  written  as  i 

V*'  -  r>  .         +  bj"fj 

5 +   d  .   -   dT   =   

f    -  f~.  '  3        f*  -  f~ 
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The  linear  scale  transformation  scheme  requires  that  the 
goals  must  be  set  within  the  range  of  (  f~,  f*  ).  Since  the 
GP  does  not  include  the  concept  of  positive  and  negative 
ideal  solutions,  we  usually  normalize  each  goal  by  way  of  a 
vector  normalization  scheme. 


Numerical  Example  2.4  s  (  continued  from  Example  2.3  ) 

When  we  set  the  equal  weights  between  the  two   goals   in 
Example  2.3,  the  non- preemptive  GP  formulation  becomes 


min    Q.5d~  +  0.5d_  1 


f1:  xl  +  di  "  di  =  30° 

0.4  0.3 

f2  =   I    .     ,   *l  +   .  X2  +  d2  -  d2 


/o.42+0.32  /o.42+0.32  -/o.42+0.32 


x1    +    x2    <    400 


2X1    +    *2    ~    500 


The  solution  is  : 
(  *x.    X2  )  *    (  250,  0  ),  d~  a  280,  d~  *  50,  and  d*  =  d+  =  0 
fjCx)  =  250   (  goal  1  not  achieved  ) 
f2(x)  =  sioo  (  goal  2  not  achieved  ) 
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CHAPTER  3   Assessment  of  Goal  Programming 

Goal  Programming  (GP)  ia  one  of  the  earliest  methods  in 
multlobjective  programming.  As  an  extension  of  linear 
programming  (LP),  GP  is  advantageous  for  its  simple 
formulation,  simple  algorithm,  and  easy  usage.  People  who 
know  LP  can  understand  GP  easily.  In  the  past  15  years, 
there  have  been  a  lot  of  articles  pertaining  to  the 
application  and  theory  of  the  GP.  The  most  recent  survey  by 
Zanakis  and  Gupta  in  1985  [Zl]  has  classified  240  GP 
articles  according  to  the  techniques  and  application  areas. 
Despite  this  fact,  the  formulation  of  GP  hide  some  major 
shortcomings  which  are  not  well-known  by  its  users  in 
solving  the  MODM  problems. 

The  most  recurrent  criticisms  relate  to  the  use  of 
preemptive  priorities  and  the  a  priori  definition  of  goals. 
GP  requires  the  decision  maker  to  Indicate  the  initial  goals 
for  each  objective  function  before  solving  the  MODM 
problems.  The  criticism  concerns  the  difficulty  of 
specifying  a  set  of  goals  a  priori  without  more  knowledge  of 
the  feasible  region.  With  regard  to  the  preemptive 
priorities,  it  is  contended  that  the  use  of  preemptive 
priorities  is  inconsistent  with  utility  function  over  the  k 
objectives,  which  ultimately  leads  to  decisions  which  are 
undesirable  to  the  decision  maker. 

The  purpose  of  this  chapter  is  an  attempt  to  address 
these  major  shortcomings  of  GP,  as   pointed   out   by   Hannan 
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[H2,H3],  Zeleny  [Z2.Z3]  and  others,  and  assess  their 
severity.  In  the  following  sections,  these  shortcomings  are 
classified,  discussed,  and  illustrated  by  numerical 
examples. 
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3.1  The  Difficulty  of  Setting  the  Initial  Goals 

A  major  problem  in  the  prior  definition  of  goals   in   GP 

is   the   difficulty   of   choosing    the    goals    as    being 

representative  of  the  decision  maker's  true  aspiration.    As 

Zeleny  [22]  mentioned,  there  are  many  financial   ratios  and 

V 
fractions  that  have  predetermined  goals   or   targets.     For 

problems   of   this   nature,   GP   is   certainly   appropriate. 

However,  in  most  practical  problems  the   value   a 

goal  should  have  is  not   apparent,   especially   in   planning 

problems,  large  scale  problems  and  public   sector   problems. 

These  problems   are   characterized   by  a   complex  decision 

process  in  which  the  appropriate  initial   goals   are   almost 

impossible  to  be  decided  by  decision  makers. 

Zeleny  [Z3]  pointed  out  that  in  most  real  life 
situations,  the  predetermined  targets,  i.e.,  initial  goals, 
are  precisely  what  decision  makers  wanted  to  learn  from  the 
analysis.  In  other  words,  he  thought  that  the  desirable 
levels  of  targets  or  goals  should  be  the  outputs,  rather 
than  the  Inputs,  of  the  problem  analysis.  Similarly,  Cohon 
and  Marks  [C5]  and  Cohon  [C6]crlticlzed  the  inconvenience  of 
specifying  target  values  for  objectives  in  GP. 

In  formulating  the  GP  model,  the  decision  maker  can  not 
just  put  the  estimated  or  expected  targets  as  the  initial 
goals  without  considering  the  constraints.  These 

estimated  or  expected  targets  are   tend   to   be   either   too 
ambitious   or   too   conservative    with    respect    to    the 


constraints.  The  GP  employs  a  minimum-distance  notion  of 
best;  that  is,  to  minimize  the  deviations  of  objective 
function  values  from  their  goals.  In  thi3  approach,  the 
optimal  GP  solution  is  quite  sensitive,  especially  the 
preemptive  GP,  to  the  goal  vector  given  by  the  decision 
maker.  If  the  goals  are  set  too  conservative,  the  solution 
space  of  GP  will  be  unreasonably  restricted  by  the  decision 
makers'  arbitrariness  and  may  result  in  a  dominated 
solution.  If  the  goals  are  set  too  ambitious,  the  unwanted 
more  weights  are  put  on  those  objectives  and  the  other 
objectives  may  suffer.  This  13  more  significant  in 
preemptive  GP.   In  most  MODM   problems,   the   objective   are 

usually  conflicting  with  each  other  competing   for   the 

limited  resources.  The  inappropriate  initial  goals  may 
result  in  the  underutilization  or  inefficient  allocation  of 
the  limited  resources. 

Although  the  decision  maker  may  have  some  past 
experiences  or  data  that  can  be  used  as  reference,  it  is 
still  very  difficult  to  decide  a  set  of  apprbpriate  Initial 
goals.  The  reason  is  that  once  the  conditions  of  problems 
are  changed,  the  coefficients,  constraints,  and  goal  levels 
of  the  model  must  be  modified.  Even  In  a  simple  GP  model  (  3 
or  4  decision  variables  ),  the  decision  maker  still  feels  it 
difficult  to  figure  out  what  variations  will  happen  to  the 
solution  space.  Unfortunately,  the  GP  approach  relies 
heavily  on  decision  maker's  perceptions  of  the  range  of 
choice  and  feasibility.   Therefore,  the  prior   determination 


of  goals  or  target  levels  without  first  exploring  the 
solution  space  could  be  too  difficult  or  too  arbitrary. 
These  inappropriate  initial  goals  will  result  in  the 
unacceptable  GP  optimal  solutions. 

The  following  numerical  example  is  a  large  scale  GP 
problem  which  was  published  by  Arthur  and  Lawrence  [A2]. 
This  paper  developed  and  solved  an  actual  problem  faced  by  a 
major  chemical  firm.  This  example  is  used  to  show  that  it 
is  not  easy  to  define  the  appropriate  initial  goals  even 
when  the  decision  maker  has  a  lot  of  past  experience  or  data 
for  reference  in  the  real  cases. 


numerical  Example  3.1  :(  A  Multiple  Goal  Blending  Problem  ) 

A.  Problem  Description  : 

The  blending  problem  is  a  variant  of  the  traditional 
mixture  problem.  This  model  is  to  decide  how  to  blend 
multiple  chemical  resources  in  order  to  produce  a  set  of 
essential  ingredients  in  such  a  manner  that  the  final 
product  will  contain  specific  percentages  of  these  essential 
ingredients.  This  problem  considers  the  blending  of  5 
raw  materials  into  10  different  final  products.  For  each 
product  and  raw  material  3  kinds  of  property  are  concerned. 
The  goals  according  to  their  priority  levels  are  listed 
below  : 
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Gl  :  The  minimization  of  the  total   coat   of   all   materials 

used  in  blending  process. 
G2  s  The  minimization  of  the  amount  of  imported  raw  material 

used  in  the  blend.  (  i.e.,  material  1  and  3  ) 
G3  :  The  minimization  of      property  2  of  all  products. 

The  minimization  of  the  impurities  in  all  products. 

The  minimization  of  the  use  of  material  4  that  requires 

special  handling  in  the  product  blend. 


G4 
G5 


B.  Model  Formulation  s 

For  the  convenience  of  discussion,  the  problem  is 
formulated  directly  from  the  data  given  in  the  paper.  This 
model  includes  50  decision  variables,  95  model  constraints, 
and  24  goal  functions  which  are  ranked  in  5  priority 
levels.  The  decision  variables  for  the  model  of  the 
blending  problem  are  : 

xij  =  tne  num*>er  of  units  of  material  1  used  in  the  blend 
of  product  j.     (  i  =   i S»  j  =  1 io  ) 


(a)  Objective  Functions  : 

Goal  1  (  Priority  1  )  :  Minimize  the  total  cost  of  all 
materials  used  in  blending  all  of  the  products.  The  initial 
goal  for  this  objective  is  9  10,000,000. 
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»4C*)  =  10  (  «X1+"ia+1'13+"14**«*,,li+"l7+"x8+"l9+*10 

*  15  (  X21+X22+X23+X24+X25+X26+X27+X28+X29+X20 

♦  5  (  *31+x32+x33+x34+x35+x36+x37+x38+x39+x30 

+  2°  (  X41+X42+X43+X44+X45+X46+X47+X48+X49+X40 

+  35  (  *51+*52+*53+x54+x55+x56+x57+x58+x59+x50 


Goal  2  (  Priority  2  )  :  Minimize  the  total  amount  of 
imported  material  1  and  3  .  The  value  in  the  parenthesl3  la 
the   initial  goal  for  the  objective. 

f2(x)  =  «11+*12+*13+*14+Xls+*i6+X17+Xt8+X19+xi0     (90,000) 
f3(x)  "  X31+X32+X33+X34+X35+X36+X37+X38+X39+X30    ("0.000) 

Goal  3  (  priority  3  )  :  Minimize  the  property  2  of  all  the 
products. 

f„(x)  =  0.02X. .+0.04x„ .+0.03x,,+0.02x. ,+0.06xc,         (2,000) 

4  —  11       21       31       41       51 

fc(x)  =  0.06x,„+0.02xo„+0.11x,„+0.05x.,+0.08xc,        (1,900) 

5  —  12        Z2        3Z        4Z        5Z 

f,(x)  =  0.032X, ,+0.024xo,+0.024x,,+0.072x. ,+0.088x_,   (1,750) 

6  —  13         Z3         33         43         53 

f-(x)  =  0.056X   +0.016X   +0.008X   +0.008X   +0.032X     (2,000) 

f„(x)  =  0.088x,I.  +  0.016x.c  +  0.032x,_  +  0.008x„c  +  0.064x  =  c.  (2,100) 
8  —  15        25        35        45        55 

f-(x)  =  0.010X,  ,+0.120x.,+0.070x,,+0.130x„,+0.060x..,  (1,475) 
y  lo  Zo  Jo  fl  6  jo 

f,„(x)  =  0.064x,,+0.008x,,+0.016x,,+0.128x.,+0.096x_,  (1,750) 
f11(X)  «  0.136x18+0.032x28+0.008x38+0.024x48+0.048x58  (1,675) 
f12(x)  =  0.060xlg+0.120x29+0.060x39+0.110x49+0.020x59   (1,950) 
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f13(x)  =  0.056x10+0.096X20+0.080x30+0.112x40+0.024x50 

(2,400) 
Goal  i  (  Priority  4  )   :   Minimize   the   impurities   in   all 
products . 

f14<*>  a  0-°3x11+0.02x21t0.05x31+0.06x41+0.04x51  (1,750) 

f15(i)  =  0-03x12+0.02x22+0.05x32+0.06x42+0.04x52  (1,750) 

f16(x)  =  0.03x13+0.02x23+0.05x33+0.06x43+0.04x53  (2,000) 

f17(x)  =  0.03x14+0.02x24+0.05x34+0.06x44+0.04x54  (2.250) 

fla(x)  =  0.03x15+0.02x25+0.05x35+0.06x45+0.04x55  (2,500) 

f19Ci)  =  0-03x16+0.02x26+0.05x36+0.06x46+0.04x  (2,000) 

f20(X)  =  0-03xi7+0-02jt27+0-05x37+0-06x47+0-04x57  (1,850) 

f21(x)  =  0.03xla+0.02x28+0.05x3a+0.06x48+0.04x58  (2,000) 

f22(x)  =  0.03x19+0.02x29+0.05x39+0.06x49+0.04x59  (2,100) 

f23(x)  =  0-03x10+6.02x20+0.05x30+0.06x40+0.04x50  (2.750) 

Goal  5  (  Priority  5  )  :   Minimize   the   use   of   material   4 
(volatile  material)  in  all  products. 

f24(^'  =  X41+X42+X43+X44+X45+X46+X47+X48+X49+X40    (8-°°°0) 
(i>)  Model  Constraints  : 
Constraint  Group  1  :  Blend  Restriction  Constraints 
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-0.030x11-0.020x21-0.000x31-0.030x41+0.060x    <  0 

(property  1  constraint  in  product  1) 

-0.040x11-0.020x21-0.030x31-0.040x41+0.000x    <  0 

(property  2  constraint  in  product  1) 

-0.0i0x11+0.020x21-0.040x31-0.030x41+0.000x51  >  0 

(property  3  constraint  in  product  1) 

-0.080x12-0.060x22-0.040x32-0.060x42-0.000x52  <  0 

-0.050x12-0.090x22-0.000x32-0.060x42-0.030x52  <  0 

-0.080x12-0.060x22+0.020x32+0.110x42-0.060x52  >  0 

-0.048x13-0.040x23-0.0i6x33-0.000x43+0.104x53  <  0 

-0.028x13-0.036x23-0.036x33-0.0i2x43+0.028x    <  0 

-0.032x13-0.032x23-0.032x33+0.054x43-0.008x53  >  0 

-0.004x14+0.012x24+0.068x34+0.012x44+0.140x54  <  0 

-0.0i6x14-0.056x24-0.064x34-0.064x44-0.040x54  <  0 

-0.004x14+0.044x24+0.020x34+0.104x44+0.152x54  >  0 

-0.032x15+0.096x25+0.064x35-0.064x45-0.056x55  <  0 

+0.076x15+0.004x25+0.020x35-0.004x45+0.052x55  <  0 

+0.124x15+0.052x25+0.040x35-0.004x45-0.136x55  >  Q 

+0.000X   +0.060X   +0.000x,,+0.190x„,+0.000x    <  0 
**■         ■i°         Jo         46         56 

-0.050x16+0.060x26  +  0.010x36+0.070x46+0.000x5  <  0 
+0.040x16+0.000x26+0.070x36+0.000x46+0.040x56  >  0 
-0.020x17+0.068x27-0.004x37+0.004x47-0.020x    <  0 
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+0.040x17-0.016x27-0.008x37+0.104x47+0.072x57  <  0 
+0.004X  +0.140X ,7+0.008x,.7+0. 092x„ ,+0. 104x  >  0 
-0.052x18-0.012x2a-0.004x3a-0.052x4a+0.020x5a  <  0 
+0.124xla+0.020x2a-0.004x38+0.0i2x48+0.036x5a  <  0 
+0.136xla+0.088x2a+0.064x38+0.136x4a+0.052x58  >  0 
-0.040x19+0.000x29+0.060x3g+0.130x49-0.070x59  <  0 
+0.020x19+0.080x29+0.020x39+0.070x49-0.020x59  <  0 
+0.070xlg+0.12Qx29+0.030x39+0.120x49+0.020x59  >  0 
+O.0O9x10+O.033x20+O.O97x30+O.0a9x4o-0.031x5o  <  0 
-0.088x10-0.048x20-0.064x30-0.032x40-0.120x50  <  0 
+0.026x10+0.098x20+0.062x30-0.046x40-0.046x50  >  0 


Conatraint  aroup.  2  :  Maximum  amount  of  each  material  used  in 
each  product  : 

x11  <  75,000 
x12  <  60,000 
x13  <  40,000 
x14  <  35.000 
x15  <  80,000 
x16  <  40,000 


x17  <  50,000 
xla  <  70,000 
xlg  <  40,000 
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x1Q  <  25.000 
x21  <  30,000 
x22  <  25.000 
x23  <  20,000 
x24  <  30.000 
x2g  <  40,000 
x2fi  <  60.000 
x2?  <  25.000 
x28  <  20,000 
x2g  <  30,000 
X2Q  <  20.000 
x31  <  30,000 
x32  <  60.000 
x33  <  50.000 
x3   <  20.000 
x35  <  40,000 
x3fi  <  75,000 
x37  <  60.000 
x3a  <  40.000 
x39  <  35,000 
X3Q  <  60,000 
x41  <  40,000 
x42  <  60,000 
x43  <  50,000 
x44  <  30,000 
x45  <  50.000 
x46  <  30.000 
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x4?  <  25,000 
x4a  <  60,000 
x49  <  40.000 
X4Q  <  70,000 
3tgi  <  50,000 
x52  <  70,000 
x53  <  60,000 
x54  <  70,000 
x55  <  50,000 
x56  <  45,000 
x5?  <  25,000 
x58  <  50,000 
x5g  <  55,000 
x5Q  <  60,000 


Conatralnt  Group   3   :   The   maximum   supply   constraint   of 
materials  : 

x1;L+x12+x13+x14+x15+x16+x17+xia+xlg+x10  <  620,000 
x21+x22+x23+x24+x25+x26+x27+x28+x29+x20  <  371,000 

*31+X32+X33+X34+X35+X36+X37+X38+X39+X30  *  611.000 
X41+X42+X43+X44+X45+X46+X47+X48+X49+X40  S  748,000 
X51+X52+X53+X54+X55+X56+X57+X5S+X59+X50  "  521'0°0 
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Constraint  Group  4  :  To  satisfy  the  minimum  demands  s 

Xll+X21+X31+X41+X51  "  40'000 
X12+X22+X32+X42+X52  ~  40'000 
X13+X23+X33+X43+X53  "  50'000 
X14+X24+X34+X44+X54  ~  45'000 
X15+X25+X35+X45+X55  ~  60'000 
X16+X26+X36+X46+X56  ~  50'000 
X17+X27+X37+X47+X57  ~  48'000 
X18+X28+X38+X48+X58  "  50'000 
X19+X29+X39+X49+X59  "  52'000 
X10+X20+X30+X40+X50  -  66'000 

In  this  problem,  the  decision  maker  set  the  initial 
goals  for  the  above  24  objectives  and  solved  the  model  by  GP 
method.  To  examine  if  these  initial  goals  are  reasonable, 
we  first  construct  the  PIS-  and  NIS-  Payoff  Table  (  Positive 
Idea  Solution  and  Negative  Idea  Solution  ) . 
Construction  of  PIS-  and  NIS-Pavoff  Table  : 
The  MODM  Formulation  : 


in   [  fl(x).  f2(x) fk(x)  J 


Max/Mir 

a-t.      x   e  X 
where   f..(x)  :  the  benefit  objective  function,  j  ■  J 
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(3.1) 


f^fx)  :  the  cost  objective  function,  i  e  I 


PIS-Pavoff  Table  : 

The  problem  is  to  find  x,  so  as  to 


and  /  or 


max   f  (x)  ,   j  e  J  (3.2) 

s.t.     x  m   x 


min  f^x)  ,   i  ■  I  (3.3) 

s.t.      x  ■  X 


NIS-Pavoff  Table  : 

The  problem  is  to  find  x,  so  as  to 

min  fj(x)  .   j  e  J                      (3.4) 

s.t.  x  e  X 
and  /  or 

max  fx(x)  ,   i  ■  I                     (3.5) 

s.t.  x  ■  X 

In  this  problem,  the  PIS-payoff  table  is  obtained  by 
(3.3)  and  the  NIS-payoff  table  is  obtained  by  (3.5).  The 
PIS.  NIS  and  the  initial  goals  are  summarized  in  Table  3.1. 
The  PIS  Is  the  best  achievable  values  for  each  objective  and 
the  NIS  is  the  possible  worst  values  for  each  objective. 
Among  the  24  Initial  goals,  we  notice  that  eight  initial 
goals  are  out  of  the  PIS  and  NIS  range.   These  Initial  goals 
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are  unreasonable.    For  example-  the  objective  f  ,  Its  worst 

value  (NIS)  can  keep  the  amount  of  property  2  of  product  5 
within  1023.31  pounds, but  the   initial   goal   set   for   this 

objective  is  2100  pounds  just   too   arbitrary   !     The 

initial  goals  for   objective   f     and   f     have   the   same 

situation.   For  other  objectives  like  f,  ,  f,,  ,  f,„  ,  f„„  , 

6     13     18     20 

and  f-,,  the  Initial  goals  were       to   ambitious.     This 

would  affect  the  achievement  of  the  other  goals  ,  especially 
when  preemptive  priorities  are  used  and  these  objectives 
have  higher  priorities. 

Another  comment  about  the  initial  goals  concerns 
the  expected  achievement  rates  of  objectives.  In  the  last 
column  of  Table  3.1  we  calculate  the  percentages  of  expected 
achievement  of  the  initial  goals  according  to  the  PIS  and 
NIS  and  find  that  the  goals  of  higher  priority  are  not  set 
at  the  higher  percentages  of  achievement.  For  instance,  the 
expected  achievement  rate  of  the  first  priority  objective 
(f1)  is  less  than  the  achievement  rate  of  the  3econd 
priority  objectives  (t^,    f  ) .   This  is   really   against   our 

intuitive  thinking.  Generally,  if  all  the  objectives  have 
the  same  weights,  the  amount  of  sacrifice  by  each  objective 
must  be  equal.  Therefore.  if  the  objective  is  really 
absolutely  more  important  than  the  other  objectives  as  in 
this  example.  Its  goal  should  be  set  at  a  higher  achievement 
rate. 
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Table  3.1  Initial  Goal  Analysis  of  Blending  Problem 


initial 

«6  of 

GOAL 

PIS 

NIS 

goal 

goal 

f  1 

6046956.67 

20633920.42 

10000000 

72.90 

f2 

34583.33 

317779.37 

90000 

80.43 

f3 

72115.38 

357577.10 

110000 

86.73 

f4 

1066.67 

4800.00 

2000 

75.00 

rs 

1250.00 

17900.00 

1900 

96.10 

f6 

2093.02 

8306.30 

1750 

105.52 

** 

n 

1980.00 

2146.67 

2000 

88.00 

ra 

560.00 

1023.31 

2100 

-232.39 

*  * 

f9 

1000.00 

8350.00 

1475 

93.54 

no 

746.67 

1996.80 

1750 

19.74 

rn 

560.00 

640.00 

1675 

-1293.75 

*  * 

f  12 

1040.00 

4400.00 

1950 

72.92 

f  13 

3002.00 

5741.64 

2400 

121.97 

*  * 

f  14 

900.00 

4000.00 

1750 

72.58 

fl5 

1376.32 

10750.00 

1750 

96.01 

f  16 

1858.14 

7217.67 

2000 

97.35 

f  17 

1237.50 

1750.00 

2250 

-97.56 

*  * 

fl8 

2624.46 

3807.69 

2500 

110.52 

»  * 

fl9 

1600.00 

6750.00 

2000 

92.23 

f20 

2086.15 

3616.00 

1850 

115.44 

*  * 

f21 

2400.00 

2800.00 

2000 

200.00 

*» 

f22 

1820.00 

4550.00 

2100 

89.74 

f23 

2048.50 

4117.54 

2750 

66.10 

f24 

55750.33 

233109.00 

80000 

86.33 

f6-fl3.fl8-f21.f22 
f8.fll-f  17 


too  ambitious 
too  conservative 
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We  conclude  that  the  ability  of  the  decision  maker  to 
set  the  appropriate  initial  goals  seems  questionable  for 
goal  problems.  Inappropriate  setting  of  initial  goals  will 
affect  the  final  solutions  a  lot.  Therefore,  the  difficulty 
of  setting  initial  goals  in  GP  will  prevent  us  to  get  the 
preferred  solutions. 
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3.2  The  possibility  of  getting  the  doainated  solutions  : 

The  second  shortcoming  of  the  GP  formulation  is  that  the 
resulting  solution  is  not  necessarily  the  "best"  one 
available  to  the  decision  maker  in  the  sense  that  it  is 
dominated  by  another  feasible  solution.  This  suboptimlzing 
feature  of  GP  is  Implied  by  the  fact  that  the  goals  are  set 
a  priori,  as  discussed  in  Hannan  [H2],  Cohon  and  Marks  [C5], 
and  Zeleny  [Z3].   This  can  be  demonstrated  by   Figure  3.1. 

In  Figure  3.1  there  are  three   objective   functions  f   , 

f2*  and  r"3  with  their  minimum  satisfactory  goal  values  g  . 
g^.      and  g^      defined    on    a    feasible    set    X.       In 

preemptive-weights  fashion,  one  might  be  interested  in 
minimizing  d±  first.  d£  second,  and  d~  last,  subject  to  the 
following  set  of  constrants  : 


^(JT)  -  tfj  +  d~ 


f2(Jf)       -  d*  +  d~2 


f3<*)  -  d+  +  rf" 


Slnce  we  want  to  reach  the  goals  as  closely  as  possible, 
we  minimize  only  the  slack  variables  d~.    d '.    and  d~ .  The 

implication  is  that  because  surplus  variables  d+,  d+,  and  d+ 
are  not  required  to  attain  their  minimum  values,  the   values 
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higher  than  g   ,    g ,    and  g     would  be  even  more   desirable.   ( 

If  we  do  not  want  to  exceed  the  goals,  we  should  also 
minimize  d   ,  d   .    and  d   .       )   At   point   A.  g        and  g        are 

satisfied  fully;  d~,  =  <C  =  °'  and  83  l3  approached  as 
closely  as  possible,  i.e.,  d~  is  minimized  on  X  subject  to 
preserving  the  zero  values  of  d~   and   d~.     Then,   point   B 


XT1 


Figure  3.1  A  priori  setting  of  goals  leads  tc 
s  u  b  0  p  t  i  m  i  z  a  t  i  0  n 


54 


represents  the  solution  to  the  preemptive  GP  approach. 
Observe  that  point  B  is  dominated  by  all  the  points  of    the 

region   BEFKB  an   embarrassing  result  !     Also,   it   is 

worth  noting  that  all  the  points  of  the  region  BEFHGCB  are 
feasible  solutions  to  this  problem  and  that  point  C 
instead  of  point  B.  may  represent  the  GP  optimal  solution  to 
the  problem  depending  on  the  coefficients  of  the  objective 
functions.  Point  C  is  also  dominated  by  all  the  points  of 
region  CIFHGC. 

Dominated  solutions  still  will  come  out  even  we   do   not 
use  the  preemptive  priorities.   Instead  of  minimizing  P.d7  + 

pzdz  +  P3<*3  ,  with  Px,  P2.  and  P3  indicating  preemptive 
priorities,  we  would  consider  the  weights  to  be  of  the  same 
order  of  magnitude  and  minimize  a  linear  combination  of 
deviations,  such  as  ^d^  +  X2d2  +  X3d3  '  tneae  ^'3  are 
archimedean  weights.  For  simplicity,  let  us  assume  X  =  \ 
=  X3-   Then  minimizing  d     +  d~   +  d~    leads  again  to  point  B  ( 

or  C  ),  a  dominated  solution.  Even  talcing  an  analog  to  the 
Loo  metric  (  MINIMAX  )  and  minimizing  the  maximum  of  the 
deviations  d^,  d^,  and  d  would  still  lead  to  a  dominated, 
suboptlmal  solution. 

After  examing  all  the  situations  above,  we  have  to 
conclude  that  the  GP  approach,  using  preemptive  or  any  other 
types  of  weighting  structure,  can  lead  to  dominated, 
suboptlmal  solutions.  The   reason   13   that   the   GP   relies 
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heavily  on  the  "principle  of  satisf icing" .  Here,  we  face  a 
classical  example  of  confusing  inputs  with  outputs  :  Goals 
should  be  the  outputs  of  the  analysis,  not  its  inputs. 
Determining  the  goals  a  priori,  without  being  able  to 
explore  the  merits  and  the  possibilities  of  the  feasible  set 
X  first,  is  one  of  the  main  reasons  for  substandard  decision 
making. 

Below  is  a  two-variable,  two  objective  numerical  example 
which  is  used  to  graphically  demonstrate  the  possibility  of 
dominated  solutions  when  GP  approach  is  used. 


I«u»erlcal  Example  3.2  :  (  Hardee  Scheduling  problem  ) 

The  Hardee  Company  makes  two  kinds  of  dolls.  Doll  A  is 
a  high  quality  one  and  Doll  B  is  of  lower  quality.  The 
respective  profits  are  30.4  and  90.3  per  doll.  Each  Doll  A 
requires  twice  as  much  time  as  a  Doll  B.  and  if  all  dolls 
were  of  type  B,  the  company  could  make  500  per  day.  The 
supply  of  material  is  sufficient  for  only  400  dolls  per  day 
(  both  A  and  B  combined  ) .  The  problem  assumes  that  all  the 
dolls  for  type  A  and  type  B  the  factory  can  make  could  be 
sold,  and  that  the  best  customer  of  the  company  wishes  to 
have  as  many  as  possible  of  type  A  doll.  The  manager 
realizes  that  two  objectives  :  (1)  the  maximization  of 
profit  and  (2)  the  maximum  production  of  doll  A,  should  be 
considered  in  scheduling  the  production. 
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Thi3  MODM  problem  can  be  formulated  as  follows  : 

(profit) 


max  f.(x)  =  0.4x   +  0.3x 


max  f2(x)  - 

xl 

subject  to 

S^C*)  - 

Xl 

+ 

X2 

< 

400 

( product  1 ) 

(material  constraint) 
(labor  constraint) 


g2(x)  =  2xx  +  x2  <  500 


x±.    x2   >   0 


where  x  and  x  are  the  number  of  Doll  A  and  that  of  Doll  B 
produced  respectively.  And  let  constraints  be  denoted  by  x 
e  X. 

The  PIS-  and  NIS-Payoff  tables  are  obtained  by  the 
formulation  of  (3.2)  and  (3.4)  to  compare  with  the  initial 
goals  in  this  problem. 


The  PIS-payoff  table  is  : 


130 
100 


100 
250* 


100 
250 


300 

0 


In  Figures  3.2  and  3.3.  point  A  is   the   ideal   solution 
for  f,(x)  and  point  E  is  the  ideal  solution  for  f,(x). 
The  NIS-payoff  table  is  z 


fl         f2 

Xl         X2 

f      o"       0 

f2     0~120        0~ 

0          0 
0        0~400 

In  Figures  3.2  and  3.3,  the  NIS  is  point  N. 

The  GP  Approach  s 

By  using  the  GP  approach,  without  first  exploring  the 
feasible  set  X,  the  manager  establishes  the  goals  for  the 
two  objectives  as  follows  : 

Gl  :  achieve  a  profit  of  at  least  990. 

G2  s  produce  at  least  180  Doll  A. 
The  GP  formulation  for  this  problem  is  : 


mi 
s.t 


n   a  =  f  (d~),  (d~)  ] 


0.4x1  +  0.3x   +  d   -  d^  =  90 


*1  +  d2  "  d2  =  18° 


x   +  x   <  400 


2x   +  x   <  500 


V  *2  *  ° 


5a 


The  solution  of  this  GP  model  is 


x1  =  180.    x2  =  60,  and   d1  =  d*  =  d~  =  dl"  =  0 


Although  both  of  the  goals  are  achieved,  the  solution  is 
dominated.  In  Figures  3.2  and  3.3,  point  G  is  the  GP 
solution  which  is  dominated  by  all  the  points  of  region 
GCEF.  After  switching  the  priority  orders  of  these  two 
objectives,  we  still  get  the  dominated  solution,  point  G. 

Instead  of  the  preemptive  priorities,  we  minimize  the 
linear  combination  of  d   and  d   and  solve  the   nonpreemptive 

problem  by  using  several  different  sets  of  weights   (X  ,X_). 

Because  of  the  different  scales  among  objective  functions, 
the  goal-norming  is  done  first,  as  stated  in  section  2.3.4. 
The  nonpreemptive  GP  model  after  normalization  is  : 


min  X1d1  +  X2d2 


S.t.      X  6  X 


0-4      „   +      0-3  -     +  .      90 

Xl  +   /    „        *2  +  dl  "  dl 


/o.32+0.42  /o.32+0.42  /o.32+0.42 


Xi    +    d2    "    d2    "    18° 


The  different  weighting  vectors,  (  0.3,  0.7  ),  (  0.5,  0.5  ), 
and  (  0.7,  0.3  )  are  tested  and  all  lead   to   the   dominated 
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solution,  x  =    225  and  x   =  0,  i.e.,  point  F  in  Figures   3.2 

and  3.3. 

MINIMAX  is  the  model  which  minimizes  the  maximum  of   the 
deviations,  slack  or  surplus,  i.e., 

min      a 

s.t.  x  m  x 

0.4                                       0.3  .-  .+  90 

x.    +    x„    +    d.    -    d.     - 


/o.32+0.42  /o.32+0.42  -/o.32  +  0.42 

Xl  +  d2  '  d2  =  18° 
d~  <  a 

dj<a 

The  solution  of  MINIMAX  still   leads   to   the   dominated 
solution,  point  G  in  Figures  3.2  and  3.3. 
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3.3  Incompatibility  with  Utility  Preferences 

This  issue  was  first  raised  by  Harrald  et  al.  [H4]  in 
developing  a  multi-goal  approach  to  resource  allocation  for 
a  marine  environmental  protection  program  for  the  U.S.  Coast 
Guard.  It  concerns  the  issue  of  whether  preemptive  GP 
is  sufficiently  flexible  to  represent  decision  maker's 
values.  On  the  contract.  Harrald  had  identified  12 
activities,  ranging  from  monitoring  the  transfer  operations 
of  tank3hlps  and  barges  to  pursuing  a  public  education 
program  that  were  surrogate  objectives  for  the  more  general 
task  of  protecting  the  marine  environment  from  discharges  of 
oil  and  other  hazardous  materials.  The  problem  was  to 
determine  the  most  effective  levels  for  these  activities 
subject  to  resource  constraints  on  personnel,  vehicles,  and 
boats. 

The  problem  was  formulated  as  a  GP  problem  with 
preemptive  weights  and  the  solution  was  not  acceptable  to 
the  decision  makers.  After  reformulating  the  problem  and 
again  arriving  at  unacceptable  solutions,  Harrald 
evaluated  the  approach.  They  concluded  that  t  -  ...  for  a 
multiple  objective  program  with  m  objectives.  having  a 
utility  function  over  the  ra  objectives  is  theoretically 
inconsistent  with  the  preemptive  GP  approach  of  specifying  a 
goal  for  each  of  the  m  objectives  and  partitioning  the  goals 
into  different  priority  classes."  This  means  that  the 
preemptive   approach   is    incompatible    with    preferences 
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representable  as  a  utility  function  and  that,  in  particular, 
changes  in  priorities  can  be  rationalized  as  a  futile 
attempt  to  state  preferences  in  a  priority  fashion. 

Sherall   and   Soyster   [S3]   formally   established    the 
connections  between  the   preemptive   and   the   nonpreemptive 
methods.   They  demonstrated  in  the  linear  case  that,  if   the 
preemptive  problem  has  an  optimal  solution,  then  there  exist- 
a  set  of  weights  for  the  nonpreemptive   problem.   such   that 
the  optimal  solution  to  the  nonpreemptive  problem  is  optimal 
to  the  preemptive   problem.      Conversely.    and   more 
Importantly,  any  optimal  solution  to  the  preemptive   problem 
is   optimal   to   the   nonpreemptive   problem.     Thus.    the 
preemptive  problem   is   subsumed   within   the   nonpreemptive 
problem  in  these  cases.   Therefore,  every  GP   problem  with 
preemptive  weights  can  be   modeled   as   a   nonpreemptive   GP 
problem,  but  the  converse  is  not  true.   There  can  be  set3  of 
weights  such  that  the   corresponding   solution   can   not   be 
obtained  by  preemptive   formulation   with  any   ordering   of 
priorities.    Decision  makers  then   run   the   risk   of   not 
obtaining  the  preferred  solution. 

Below  is  a  simple  numerical  example  used  to  construct 
the  set  of  equivalent  weights  and  demonstrate  the 
incompatibility  between  preemptive  and  nonpreemptive 
approach. 

numerical  Example  3.4  :  (modified  Hardee's  scheduling  problem) 
The  Hardee  Company  makes  two  kinds  of  dolls.   Doll  A   is 
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a  high  quality  one  and  Doll  B  i3  o£  lower  quality.  The 
respective  profits  are  90.4  and  90.6  per  doll.  Each  Doll  A 
requires  twice  as  much  time  as  a  Doll  B.  and  if  all  dolls 
were  of  type  B,  the  company  could  make  500  per  day.  The 
supply  of  material  is  sufficient  for  only  400  dolls  per  day 
(  both  A  and  B  combined  ).  The  problem  assumes  that  all  the 
dolls  for  type  A  and  type  B  the  factory  can  make  could  be 
sold,  and  that  the  best  customer  of  the  company  wishes  to 
have  as  many  as  possible  of  type  A  doll.  The  manager 
realizes  that  two  objectives  :  (l)  the  maximization  of 
profit  and  (2)  the  maximum  production  of  Doll  A,  should  be 
considered  in  the  scheduling  the  production. 

By  using  the  GP  approach,  the   manager   establishes   the 
goal3  for  the  two  objectives  as  follows  : 

Gl  :  achieve  a  profit  of  at  least  9250. 

G2  :  produce  at  least  200  Doll  A. 
The  constraint  and  goal  functions  are  shown  in  Figure  3.4. 

Preemptive  Approach  : 

Min  a  =  T  (d~),  (d~)  1 

Subject  to 

0.4x1  +  0.6x2  +  d~  -  d*  =  250         (3.6) 

Xl  +  d2  "  d2  '    200 
*X    +    x2    -    400 
2x1  +  x2  <  500 
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Xl'  X2'  dl'  dl'  d2'  d2  -  ° 


(1)  The  preemptive  solution  for  problem  (3.6)  is  point  A  : 

Xl  =  °'  x2  =  400'  dl  =  10'  d2  =  200'  and  dt  =  d2=  °" 

(2)  The   preemptive   solution   for   problem   (3.6)   with 
priorities  Interchanged  is  point  C  : 

xx  =  200.  x2  =  100.  d~  =  110.  and  d*  =  d~  =  d*  =  0. 


Nonpreemptlve  Approach  : 

If  \^  and  X2  represent  the  weights  of  the  two  goals  and 
assuming  that  the  goal  norming  has  been  considered  in 
assigning  the  weights.  the  nonpreemptive  GP  formulation 
would  be  : 

Min  a  =  [  ^d-  +  X2d"  ] 

Subject  to 

0.4x1  +  0.6x2  +  d~-  d*  =  250  (3.7) 

Xl  +  d2  "  d2  =  200 
X1  +  x2  <  400 

2x1    +  x2  <  500 

Xl'  V  V  *t'    d2'  d2  "  ° 


500 


400 


300 


200 


100 


Fdgure  3.4  Decision  variable  space  representation 
of  modified  Hardee  Co.  problem 


After  conducting  the  sensitivity  analysis.   we   get   the 
following  solutions  :  (  if  x   +  X   =  1  ) 

(1)  if  0  <  X1  <  5/9,  the  nonpreemptive  solution  is  point  C  : 
xx    =    ZOO.  x2  =  100,  d~  =  110,  and  d+  =  d~  =  d+  =  0. 
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(2)  if  5/9  <  X1  <  5/6,  the  nonpreemptive  solution  is  point  B 
xx  =  100.  x2  =  300,  d~  =  30,  d~  =  100,  and  d+  =  d+  =  0. 

(3)  if  5/6  <  \1  <  1,  the  nonpreemptive  solution  is  point  A 
Xl  =  °*  x2  =  400'  dl  =  10'  dz  =  200'  and  di  =  d2  =  °* 


Thus,  the  comparison  of  the  results  for   the   preemptive 
and  nonpreemptive  approach  is  i 


0   *  *•!  -  5/9     same  solution   as   preemptive   problem 
(3.6)  with  priorities  interchanged 

5/9  <   *j  *  5/6     no  preemptive  solution  with   same  set  of 
goals  will  yield  the  same  solution 

5/6  -  xi  *   x      same  solution  as   preemptive   problem 
(3.6) 


If  we  think  of  the  objective  of  the  original  problem  as 
consisting  of  a  linear  utility  function  of  the  goals  with 
unknown  weights,  then  one  preemptive  order  will  generate  the 
same  solution  as  a  range  of  weights  (  0  <  \  <  5/9  )  and  the 
other  possible  preemptive  order  will  generate  the  same 
solution  as  another  range  of  weights  (  5/6  <  X  <  1  ). 
However,  if  the  actual  utility  function  is  in  a  third 
range  (  5/9  <  ^  <  5/6  ),  then  neither  preemptive 
formulation  will  yield  the  solution.   This  also  proves   that 


66 


changing  the  priorities  of  the  goals  alone   cannot   yield   a 
good  approximation  to  a  utility  function. 

Significantly,  the  preemptive  approach  is  worse  than  the 
nonpreemptive  approach.  But  this  does  not  mean  that  the 
nonpreemptive  GP  is  good  enough  to  approximate  the  utility 
function  of  decision  maker's  preference.  One  problem  in 
using  this  approach  is  the  well-known  difficulty  in 
specifying  a  priori  the  initial  goals  and  the  deviational 
weights.  Another  problem,  not  generally  recognized,  is  the 
tendency  for  optimal  solutions  to  standard  linear  goal 
programs  to  occur  at  extreme  points.  One  would  not  usually 
expect  the  activity  levels  maximizing  a  utility  function  to 
be  at  an  extreme  point,  especially  if  the  extreme  points  are 
not  "close"  to  each  other,  such  as  the  extreme  points  in 
Figure  3.4. 
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3.4  'Absolutely  more  important"  is  not  common 

In  1965,  Ijiri  [13]  presented  a  definition  of  preemptive 
priority  levels  so  as  to  treat  goals  according  to  their 
perceived  importance  and  solve  the  problem  of 
incommensurable  objectives.  Although  GP  need  not  include 
the  preemptive  priority  concept,  the  idea  has  been  so  widely 
used  that,  when  most  people  now  talk  about  GP.  they  are 
really  referring  to  GP  within  a  preemptive  priority 
structure. 

There  are  situations  in  which  a  decision  maker  feels 
that  the  objectives  are  clearly  hierarchically  ordered.  It 
will  be  impossible  for  him  to  give  information  about 
trade-offs.  because  the  one  objective  is  so  important 
compared  with  the  other  that  they  can  not  be  trade-off 
against  one  another.  In  this  situation,  if  the  targets  are 
also  very  clear,  the  preemptive  GP  should  be  used.  But  in 
real  application,  this  kind  of  situation  is  not  common. 

In  the  survey  of  the  GP  application,  the  GP  model  is 
usually  formulated  into  many  preemptive  priority  levels  as 
compared  with  the  number  of  objectives  (  goals  ).  It  is 
just  unrealistic  in  the  practical  problems  that  the  goals 
can  be  ranked  in  so  many  preemptive  levels.  It  seems  the 
fact  is  that  people  try  to  use  the  preemptive  GP  algorithm 
as  a  surrogate  :  (1)  for  problems  with  commensurable  goals 
to  avoid  specifying  weights  for  the  goals  and  (2)  for 
problems   with   incommensurable   goals   to   avoid   different 
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scales.   As  In  Lee's  text  [LI]  : 

"...  question  that  may  arise  is  whether  the  preemptive 
priority  structure  on  an  ordinal  ranking  is  better  than 
utility  analysis.  If  there  existed  a  systematic  methodology 
that  could  effectively  measure  the  decision  maker's  utility 
for  multiple  goals,  obviously  the  utility  approach  would  be 
most  convenient..." 

"...  However,  there  is  no  effective  utility  analysis 
method  at  present.  Furthermore,  under  normal  circumstances 
the  decision  maker  is  not  capable  of  measuring  precisely  how 
much  more  important  the  first  goal  is  than  the  second  in 
cardinal  value.  Often  the  only  knowledge  the  decision  maker 
has  is  that  In  his  judgement,  a  first  goal  is  more  important 
for  the  organization  than  a  second  goal.  Hence.  a 
preemptive  priority  structure  in  terms  of  ordinal  ranking  of 
goals  is  the  more  appropriate  vehicle  for  decision 
analysis." 

It  is  true  that  people  can  take  account  of  only  a  very 
limited  number  of  criteria  in  making  an  evaluative 
judgement.  But  "more  important"  does  not  equal 
"absolutely  more  Important  --  a  lexicographic  order".  On 
the  other  hand,  if  the  utility  approach  is  too  difficult  to 
use.  it  does  not  mean  that  we  have  to  choose  the  preemptive 
approach.  There  are  still  some  other  MODM  techniques,  such 
as  the  interactive  approach, that  do  not  have  many  of  the 
practical  and  theoretical  shortcomings  of  the  noninteractive 
methods.  The  use  of  the  preemptive  approach  will   have   some 
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drawbacks  for  solving  the  problems  of  not  having  the 
preemptive  importance  in  relationship  to  one  another  among 
objectives.  Except  the  one,  discussed  in  previous  section, 
that  the  preemptive  weight  is  not  compatible  with  the 
utility  function.  Also,  the  rigidity  of  lexicographic 
orders  allows  no  trade-off  between  objectives  of  different 
priority  levels. 

Gass  [Gl]  had  the  following  comment  about  the  preemptive 
GP  approach  :  "  ...  I  have  always  felt  that  such 
formulations  hide  the  fact  that  they  greatly  reduce  the 
range  of  acceptable  compromise  solutions  by  sequentially 
cutting  off  portions  of  the  allowable  (feasible)  solution 
space  ...-.  The  restrictive  nature  of  preemptive  goals 
takes  away  the  opportunity  to  Investigate  a  reasonable  set 
of  possible  solutions  from  which  a  final  solution  is 
selected.   This  can  be  demonstrated  in  Figure  3.5. 

In  Figure  3.5  there  are  two  objective  functions   f    and 
f2  with  their  minimum  satisfactory  goal   values   g    and   g 
defined  on  a  feasible  set  X.   In  preemptive-weight  approach, 

one  might  be   Interested   in  minimizing   d~   first   and   d~ 

1  2 

second,  or  d£  first  and  d~  second,  subject  to  the   following 
set  of  constraints  s 

f2Cx)       ♦  d^         -d*    =  g2 


If  we  minimize  d  first,  the  feasible  region  NABGH  will 
be  cut  off  in  the  first  iteration  of  calculation,  and  point 
B  will  be  the  solution  to  the  preemptive  GP  approach  when  we 


minimize 


in  the  second   Iteration. 


In   the   same  way. 


Figure  3.5  The  restrictive  nature  of  preemptive 
Goal  Programming 


another  possible  solution  will  be  point  C  if  we  minimize 


first.   No  matter   if 


d   or  d   are  minimized   first 


cut  off  large  portions  of  solution  space  and  the  feasible 
solutions  on  the  line  section  BC  are  never  considered.  Even 
worse,   once   the   goal3   of   higher   priority  are  set   too 
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ambitious  and  can  not  be  satisfied,  such   as   f '   in   Figure 

3.5,  then  the  unique  solution  (  point  F  )  is  the  preferred 
solution  to  the  entire  problem, and  the  objectives  of  lower 
priority  would  not  be  considered  at  all. 

Another  danger  of  preemptive  GP  ia  that  the  preemptive 
GP  explicitly  assumes  that  goals  within  the  same  priority 
level  should  be  commensurable  (  comparable  ) .  Due  to  this 
presumption,  Rlnguest  and  Gulledge  [R2]  warned  that  goals 
which  were  of  equal  Importance  might  be  placed  in  different 
levels  only  because  they  were  incommensurable. 

Below,  there  is  a  numerical  example  from   Lee  [LI]   for 

illustrating  how  people  utilize  the  preemptive  GP  to   solve 

the  MODM  problem  and  what  the  disadvantages  are  for  this 
approach. 


Numerical  Example  3.5  :  (  Continental  Electronics.  Inc.  ) 

Continental  Electronics.  Inc.  produces  two  products  : 
record  players  and  tape  recorders.  The  production  of  both 
products  is  done  in  two  separate  machine  centers  within  the 
plant.  Each  record  player  requires  two  hours  in  machine 
center  1  and  one  hour  in  machine  center  2.  Each  tape 
recorder,  on  the  other  hand,  requires  one  hour  in  machine 
center  1  and  three  hours  in  machine  center  2.  In  addition, 
each   product   requires   some   in-process   inventory.     The 
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per-unit  in-process  Inventory  required  is  350  for  the  record 
player  and  9  30  for  the  tape  recorder. 

The  firm  has  normal  monthly  operation  hours  of  120  for 
machine  center  1  and  150  for  machine  center  2.  The  cost 
accounting  department  estimates  that  the  average  costs  of 
per -hour  operation  are  380  and  320  for  machine  centers  1  and 
2, respectively.  The  estimated  profit  per  unit  is  3100  for 
the  record  player  and  975  for  the  tape  recorder.  According 
to  the  marketing  department,  the  forecasted  sales  for  the 
record  player  and  the  tape  recorder  are  50  and  80. 
respectively,  for  the  coming  month.  These  data  can  be 
summarized  into  the  following  table  : 


Record 
Player 

Tape 
Recorder 


Machine   Machine   In-process   Profit      Forcasted 
center    center    Inventory   per  unit     sales 
*    1       #2       cost 


2  hr. 


1  hr.     3  50 


1  hr.     3  hr. 


3  30 


3  100      50  units 


3  75 


80  units 


Machine 
center  #1 


Normal  operation 
hours  /  month 


120  hr. 


Average  cost  of 
operation  /  hour 

3  80 


Machine 

center  #2 


150  hr. 


3  20 
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The  president  of  the  firm  has  established  the  following 
goals  for  production  in  the  next  month,  in  ordinal  ranking 
of  importance  : 

Gl  :  Limit  the  amount  tied  up  in  in-process  inventory  for 

the  month  to  94,600. 
G2  :  Achieve  the  sales  goal  of  50  record  players  for  the 

month . 
G3  :  Avoid  any  underutilization  of  regular  operation  hours 
of  both  machine  centers.  (  The  president  desires  to 
assign  weights  according  to  the  cost  of  idle  time  in 
each  machine  center.) 
G4  :  Limit  the  overtime  operation  of  machine  center  i  to 

20  hours. 
G5  s  Achieve  the  sales  goal  of  80  tape  recorders  for  the 

month. 
G6  :  Limit  the  sum  of  overtime  operation  for  both  machine 
centers   (the   president   wants   to   assign   weights 
according  to   the   overtime   cost   of   each  machine 
center). 
A.  Formulation  of  GP  Model  : 

The  complete  production  planning  model  can  be  formulated 
as  below  : 

Mm  z  =  [  Pid*  ,  P2a;  .  pi(4d;  +  d;)  ,  P„dt,  .  P=d: 


■2~4    '    r3l~l    T   u2;     '    r4Qll     '    r5a5 


p6(4d 


!*<>       ] 
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2X1  +    X2  +  dl  "  dl  =    12° 

Xl  +   3X2  +  d2  "  d2  *    150 
50x1  +  30x2  +  d^  -  d*  =  4,600 

xl         +  d4  "  d4  =     50 

X2         +  d5  "  d5  *     30 

dt  +  dll  "  dtl  =     20 

X1'X2  "  °' 

d~.  d±  >  0,  i  =  1,2,3,4,5.  and  11 

where  xA  =  number  of  record  players  produced  in  the  month 

x2  =  number  of  tape  recorders  produced  in  the  month. 

B.  Solutions  : 

The  solution  of   the   above   problem  of   preemptive   GP 
formulation  is  : 

x1=50,  x2=40,  d+=20,  d2=20,  d~=900,  d~=40, 

dI  =  d2  "  d3  "  d4  =  dI  =  d5  =  dtl  "  Q 
The  solution  indicates  that  the  firm  produces  50  record 
players  and  40  tape  recorders  with  20  hours  of  overtime 
operation  in  both  machine  centers.  The  total  in-process 
inventory  is  $3,700,  3900  short  of  the  94,600  limit  as  set 
by  the  president.  The  firm  is  not  able  to  achieve  the  sales 
goal  of  30  tape  recorders.  The  total  profit  of  the  firm  in 
the  month  is  38,000. 
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C.  Analysis  and  Discussion  : 

In  this  problem,  there  are  6  preemptive  priority  levels. 
Although  the  flr3t  priority  level  goal3  are  usually  looked 
as  the  physical  constraints  in  preeemptive  GP  model,  it  i3 
still  not  realistic  that  the  goals  can  be  ranked  in  so  many 
preemptive  priority  levels.  Some  goals  may  be  more 
important  than  the  other,  but  the  relative  weights  should 
not  be  all  in  preemptive. 

If  it  is  true  that  the  formulation  of  preemptive  GP  in 
this  problem  is  just  a  tool  to  solve  the  problem.  The 
trade-off  among  objectives  should  be  the  problem,  the 
slight  modification  of  thi3  problem,  the  first  3  goals  are 
considered  to  be  the  physical  constraints  and  a  tradeoff 
among  the  other  goal3  i3  compared.  Thi3  new  MODM  problem 
will  be  formulated  as  follows  : 

Max  fjCx)  =  100x1  +  75x2    (profit) 

Mln  f4(x)  =  2xx    +  x2     (G4:  overtime  of  m/c  center  1) 

Max  f5(x)  =  x2  (G5:  tape  recorder) 

s.t. 

50x1  +  30x2  <  4600    (Oil  in-proce33  inventory  cost) 


player ) 
2xx  +    x2  -   120    (G3'  no  underutilization  m/c  1) 

*1  +   3x2  ~   150    <G3:  no  underutilization  m/c  2) 

V  x2  *  ° 


76 


In  Figure  3.6.  the  feasible  region  for  this  new 
formulation  i3  BCD  and  the  PIS-  and  NIS-payoff  tables  are 
obtained  as  follows  : 

PIS-Payoff  Table  i 


10250 


170 


f5    10250 


170 


70 


f4     7500    133.33    33.33 


70 


50 
50 
50 


70 
33.33 

70 


In  Figure  3.6,  point  C  is  the  PIS  for  f,(x)   and   f„(x), 

1  —  5  — 


and  point  B  is  the  PIS  for  f„(x). 

4  — 


NIS-Pavoff  Table  : 


fl 

f4 

f5 

Xl 

X2 

fl 

7500" 

133.33 

33.33 

50 

33.33 

f4 

9562.5 

179.17" 

24.17 

77.5 

24.17 

f5 

9562.5 

179.17 

24.17" 

77.5 

24.17 

In  Figure  3.6,  point  B  is  the  NIS  for  f  (x)  and  point   D 
is  the  NIS  for  f^Cx)  and  fg(x). 


130 


100 


5  Ox. +30x=4600(inventory    cost) 
2x1+x2»120(m/c    1) 


x    =50(product    1) 


j+jXj-lSOCa/c  2) 


Figure  3.6  The  decision  variable  space  representation 
in  Continental  Electronic,  Inc.  problem 

From  the  PIS-payoff  table.  we  can  see  the  maximum 
profit  that  can  be  made  is  910,250  (  as  compared  with 
38.000  la  the  preemptive  GP  solution  ).  and  the  maximum 
amount  of  tape  recorders  is  70  units  (  as  compared  with 
40  units  In  the  preemptive  GP  solution),  both  under  the 
sacrifice  of  the  overtime  operations  of  machine  centers. 
Although  these  positive  ideal  solutions  may  not  be  the 
preferred  solution  to  the  decision  maker  for  the  high 
overtime  operations  of  machine  centers,  this  demonstrates 
that  the  decision  maker  will  have   more   choices   about   his 
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preferred  solution  if  the  preemptive  GP  is  not  used. 

In  this  example,  if  a  tradeoff  between  the  overtime 
operation  of  machine  centers  and  the  amount  of  the  tape 
recorders  can  be  made,  we  can  also  make  more  profit.  For 
instance,  the  solution  x1=  50  and  x£  =  52  will  have  the  same 

achievement  rate  of  60%  for  both  objectives  f  and  f 
according  to  the  model  and  the  profit  for  this  solution  is 
38,900,  this  is  point  F.  Even  better,  if  we  can  trade-off 
between  the  amount  of  tape  recorder  and  record  player,  we 
can  make  more  profit  without  increasing  the  overtime 
operation  of  machine  centers. 

The  basic  approach  of   dealing   with   the   MODM   problem 
should  be  a  kind  of  compromise  among  the  objectives   to   get 
the   preferred   solution.     The   preemptive   GP.    however, 
disallows  the  tradeoffs  of  small  losses  in   higher   priority 
objectives  for  larger  gains  in  low  objectives.   Even  worse, 
the   analyst   is   usually   advised   to   perform   sensitivity 
studies  to  investigate  the  impact  of   switching   priorities. 
This  also  proves  that  there  is  no  preemptive  priority  to  the 
goals   in   decision   maker's   mind.     The   reason   to    use 
preemptive  weights  is  just  because  the  preemptive  GP  is  easy 
to  formulate  and  the  algorithm  is  easy  to  understand.     The 
improper  use  of  this  weighting  scheme  can   cause   more   harm 
than  good.   Therefore,  GP  with  preemptive  weights  is   not   a 
recommended   approach   unless   the   goals   really   do    have 
preemptive  importance  in  relationship  to  one  another. 


CHAPTER  4   Enhancement  of  Goal  Programming 

In  Chapter  3,  we  discussed  many  shortcomings  of  Goal 
Programming.  But  it  does  not  mean  that  Goal  Programming 
should  not  be  used  in  any  circumstances.  Actually,  there 
are  some  situations  in  which  Goal  Programming  is  more 
suitable  than  other  MODM  methods.  Therefore,  if  the  GP  is 
chosen  as  the  desirable  solution  technique,  it  would  be  nice 
to  have  some  techniques  to  remedy  or  prevent  these 
shortcomings.  In  this  chapter,  techniques  for  this  purpose 
are  suggested.  And  a  general  decision  rule  is  also 
presented  to  help  the  user  to  identify  whether  the  GP  is  the 
right  method  to  be  used.  If  the  GP  is  not  the  appropriate 
method,  other  MODM  methods  are  recommended. 


4.1   Remedial  Methods  for  Goal  Programming 

In  this  section,  four  remedial  methods  are  suggested  and 
discussed.   These  remedial  methods  include  : 

(1)  Incorporation  of  PIS,  NIS,  MAR,  and  MAG 

(2)  Techniques  for  checking  dominated  solutions 

(3)  Use  of  non-preemptive  weights 

(4)  The  interactive  approach  of  GP 

These  methods  will  be  explained  in  detail  in   the   following 
subsections. 
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4.1.1  Incorporation  of  PIS.  MIS.  MAR  and  HAG 

In  GP.  it  is  difficult  for  decision  maker  to  state  the 
goals  without  any  knowledge  about  the  feasible  objective 
space.  Except  in  very  few  cases,  the  initial  goals  of  GP 
approach  usually  roughly  set.  Every  so  often,  the 
initial  goals  are  solved  to  be  inappropriate  or  meaningless 
with  respect  to  the  problem  formulation.  In  Chapter  3,  we 
discussed  that  the  Inappropriate  initial  goals  can  result  in 
dominated  solutions  and  inefficient  allocations  of 
resources.  So,  it  would  be  a  good  suggestion  for  GP 
approach  to  first  explore  the  feasible  objective  space  and 
give  the  decision  maker  some  idea  about  what  kind  of 
compromise  solutions  can  be  obtained . 

In  this  section,  we  introduce  the  concepts  and 
algorithms  of  the  Positive  Ideal  Solution  (PIS),  the 
Negative  Ideal  Solution  (NIS),  the  Maximum  Achievable  Rate 
(MAR),  and  the  Maximum  Achievable  Goal  (MAG).  These 

can  serve  as  excellent  reference  values  to  decision 
makers  in  setting  the  Initial  goals.  The  PIS  and 
NIS  provide  the  upper  and  lower  bounds  of  the  objective 
function  values.   The  set  of   MAG   is   a   feasible   solution 

which  has  the  highest  same  rate,  MAR,   for   all   objectives. 

These   concepts   have   been   used   in   Hwang    and    Masud's 

Interactive  Sequential  Goal  Programming  II  (ISGP   II)   [H7]. 

The  meaning  of  MAR  and  MAG  have  been  stated  in  [H7]  as  : 

" In   a   group  decision   making   problem,   if   each 

objective  represents  one  decision  maker's  payoff,   then   the 
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positive-ideal  solution,  f*,  if  attainable,  would  make  each 
decision  maker  happy  because  it  simultaneously  optimizes 
each  decision  maker's  payoff.    However.   in  a   real  world 

* 
problem-  f   is  usually  not  attainable.     To  dissolve   the 

group  conflict,  a  compromise  among  the  decision  makers  is 
necessary.  That  is.  each  decision  maker  expects  to  achieve 
only  portion  of  the  f*  instead  of  the  entire  f*.  The  amount 
of  sacrifice  by  each  decision  maker  must  be  equal  in  order 
to  maintain  a  "cooperative  group  spirit".  Hence.  MAR  and 
MAG  are  introduced  to  satisfy  that  condition.  ...    - 

The  MAG  can  be  interpreted  as  a  set  of  nondominated 
solutions  in  which  each  objective  has  the  maximum  percentage 
of  achievement  (MAR)  with  respect  to  the  feasible  objective 
space. 

The  algorithms  to  solve  the  PIS.  NIS.  MAR,  and  MAG  will 
be  described  below  : 


Problem  Formulation  : 

The  MODM  problem  Is  : 

Max/Min  [f^x),  f^*) f^x)    J  (4  ^ 

3-t.    x  e  X 

where  f j(X)  :  the  benefit  objective  for  maximization,  j  e  j 
f±(x)  :  the  cost  objective  for  minimization.  1  <=  I 
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The  Positive  Ideal  Solution  (PIS)  (  f *  )  : 

The  positive  -ideal  solution  is  the  one   that   optimizes 
each  objective  function  simultaneously,  e.g., 


max   f  (x)  ,   j  e  J 
xex    3 


(4.2) 


and/or 


min   f  (x)  ,   i  e  I 
xeX    * 


(4.3) 


An  ideal  solution  can  be  defined  as  f*  =  (  f*,  f*.  .     f*  i 

1'   2'  ■  '  •  '   fc 

where  fg  is  a  feasible   and   optimal   value   for   the 

objective  function.   The  PIS-Payoff  Table  is  given  as 


,th 


PIS-Payoff  Table 


'21 


kl 


12 


"k2 


'Ik 


"2k 


The  Negative  Ideal  Solution  (NIS)  (  f ~  )  : 

Simllary,   the   negative-ideal   solution   is   the   worst 
possible     solution     for     each     objective     function 
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simultaneously,  e.g.. 


min  f  .(x)  .   jeJ 
x«X   3 


(4.4) 


and/or 


max  f  ( x )  .   lei 
xsX 


(4.5) 


ffc  )   where 


A  NIS  can  be  defined  as  f~  =  r  f~   f" 

1   1'   2'  -  " 

f^  is  a  feasible  and  the  worst  value  for  the  _£th  objective 

function.   The  PIS-Payoff  table  is  given  as  : 


PIS-Payoff  Table 


fl     Z12 


221    f2 


kl     k2 


'Ik 


2k 


The  Maximum  Achievable  Rate  (MAR,  B%)  and  Maximum  Achievable 
Goal  (MAG,  fg)  s 

The  problem  is  to  find  x  so  as  to  : 
k 


min  a 


9  . 


[  I    "A   ] 


(4.6) 


r  =  l 


84 


3.t     x  e  X 

£.,(*)  +  (  <*~.  -  a*  )(  r*  -  f~  )  =  b*  .     j  e  j 

f±(x)  -  (  d"  -  <J+  )(  f^  -  £*  )  *  b*  .   i  e  I 

0  -  dr  '  dr  -  1'  r  =  1'2 'k 


where   b*  =  t *  -  (  1  -  B%  )(  { *   -  f"  ) ,    j  e  j 


b 


i 


f±  +  (  i  -  b%  >(  f"  -  e*  ),   i  e  ! 


The  MAR  is  unknown  initially.  To  properly  3tart  the 
calculation  of  (4.6),  we  have  to  a33ign  a  value  to  B%.  The 
usual  value  of  B%  to  start  with  is  75%.  The  following  steps 
illustrate  the  algorithm  to  find  the  MAG  and  MAR. 

Steg  3^1  :  Check  if  fg  is  a  set  of  nondorainated  solutions. 

If  a   =  0,  go  to  step  3.2. 

If  a  *  0,  we  decrease  B%  by  5%.  i.e.,  B%  =  (B  -  5)%, 
recalculate  (4.6).  and  check  the  value  of  ag  again.  The 
procedure  is  repeated  until  an  answer  a9  =  o  is  found. 
After  an  a  =  0  is  found,  we  Increase  B%  by  1%,  i.e.  B%  =  (B 
+  1)%,  recaluate  (4.6),  and  check  the  value  of  a9  again. 
The  procedure  repeats  until  an  a9  *  0  is  found,  and  that  ag 
*   0  signals  that  MAG  and  MAR  are  found. 
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Step  3.2  : 

When  a  =  0,  we  increase  B%  by  5%.  i.e.,  B\  =  (B  +  5)%, 
recalculate  (4.6),  and  check  the  value  of  ag  again.  The 
procedure  is  repeated  until  an  answer  ag  *  0  is  found. 
After  an  a  »  0  la  obtained,  we  decrease  B%  by  1%,  i.e.  B%  = 
(B  -  1)%,  recalculate  (4.6),  and  check  the  value  of  ag 
again.  The  procedure  is  repeated  until  an  answer  ag  =  o  is 
found,  and  that  ag  =  0  signals  that  MAG  and  MAR  are  found. 
We  conclude  that  MAR  is  (B  +  1)%  and  MAG  is  the  set  of 
nondominated  solutions  obtained  at  (B  +  1)%. 


The 


numerical  Example  4.1  :  (  Hardee  Scheduling  problem  ) 

The  Hardee  Company  makes  two  kinds  of  dolls.  Doll  A  i 
a  high  quality  one  and  Doll  B  is  of  lower  quality, 
respective  profits  are  90.4  and  30.3  per  doll.  Each  Doll  A 
requires  twice  as  much  time  as  a  Doll  B,  and  if  all  dolls 
were  of  type  B,  the  company  could  make  500  per  day.  The 
supply  of  material  is  sufficient  for  only  400  dolls  per  day 
(  both  A  and  B  combined  ).  The  problem  assumes  that  all  the 
dolls  for  type  A  and  type  B  the  factory  can  make  could  be 
sold,  and  that  the  best  customer  of  the   company   wishes   to 
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have  as  many  as  possible  of  type  A  doll.  The  manager 
realizes  that  two  objectives  :  ( 1 )  the  maximization  of 
profit  and  (2)  the  maximum  production  of  doll  A,  should  be 
considered  in  scheduling  production. 

The  MODM  problem  can  be  formulated  as  follows  : 

max  f  (x)  =  0.4x.  +  0 . 3x 
1  12 

max  f2(x)  =  x^ 
s.t. 


x1  +  x2  <  400 


2x1  +  x2  <  500 


V  x2  *  ° 


where  x1  and  x£  are  the  number  of  Doll  A  and  that  of 
Doll  B  produced,  respectively.  Let  constraints  be  denoted  by 
X   e  X. 


The  PIS  and  NIS 


We   solve   max   f^x)   and   max   f.,(x)   subject   to   the 
constraints,  respectively,  and  obtain  the   PIS-Payoff   Table 


as 
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fl 

f2 

fl 

* 
130 

100 

f2 

100 

* 
250 

£   =  (  130,  250  ) 
Similarly,  we  solve  min  f^Cx)  and  min  f£(x)  subject  to  x 
X  respectively,  and  obtain  the  NIS-Payoff  Table  as  : 


fl 

f]k  o"       0 

f2       o       o" 


£  =  (  o,  o  ) 

The  MAR  and  MAG  : 

We  set  the  B%  =  75%,  and  solve  the  following  problem  : 
min  a9  =  <  d^  +  d~  }  (4.7) 

s.t. 

x  m   x 


fjCx)  +  130  (d~  -  d+)  =  bj 
f2(x)  +  250  (d~  -  d  +  )  =  b2 


0  <  d  ,  d,  <  1,  r  =  1,2 
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where  b1    =    130  -  (  1  -  75%  )(  130  -  0  )  =  97.50 
b2  =  250  -  (  1  -  75%  )(  250  -  0  )  =  187.50 
The  solution  la  given  aa  : 

x  =  (  187.50.  75.00  ) 
fg  =  (  97.50,  187.50  ) 


The  answer  a9  =  0  implies  that  with  the   given   b   there 

exlst3  a  dominated  f   which  la  not  appropriate  to   serve   aa 
MAG. 


Since  B%  =  75%  results  in  a  dominated  f9,  we  increase  B% 
by  5%,  i.e.,  B%  =  80%,  and  recalculate  (4.7)  with 

b1  =  130  -  (  1  -  80%  )  (  130  -  0  )  =  104.00 

b2  =  250  -  (  1  -  80%  )  (  250  -  0  )  =  200.00 

The  solution  ia  given  aa  J 

x  =  (  200.00,  80.00  ) 

f9  =  (  104.00,  200.00  ) 

a9  =  0 
Again,  we  find  the  new  ag  =  0,   i.e..   f9   ia   a   aet   of 
dominated  aolutions.   We  can  increase  B%  by  another  5%,   i.e. 
B%  -  85%.   Another  calculation  of  (4.7)  ia  conducted  with 


b 


1 


130  -  (  1  -  85%  )  (  130  -  0  )  =  110.00 


b2  =  250  -  (  1  -  85%  )  (  250  -  0  )  =  212.50 
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The  solution  is  given  as  : 

X  =  (  212.50,  75.00  ) 
f9  =  (  107.50,  212.50  ) 
ag  =  0.0231 

The  answer  ag  *  0  implies  that  fg  is  a  nondominated  solution 
The  MAR,  B%,  is  located  somewhere  between  85%  and  81%.  To 
locate  the  precise  B%,  we  decrease  B%  by  1%.  i.e.  B%  =  84%, 
and  calculate  (4.7)  with  b1  =  109.20  and  b2  =  210.00.  The 
solution  is  gliven  as 

x  =  (210.00,  80.00) 

fg  ■  (108.00,  210.00) 

a9  =  0.0092 
Since  the  answer  ag  *  0,  we  decrease  B%  by   another   1%, 
i.e.  B%  =  83%,  and  proceed  to  another  calculation   of   (4.7) 
with  b  =  (107.90,  207.50).   The  solution  is  given  as 

x  =  (207.50,  83.00) 
f9  =  (107.90,207.50) 
ag  =  0 

The  answer  ag=  0  implies  that  the  fg,  when  B%  =  83%,  is  a 
set  of  dominated  solutions.  We  conclude  that  MAG  is  located 
at  fg  =  (108.00,  210.00)  and  MAR  is  84%,  i.e.  B%  =  84%. 

The  results  can  be  presented  to  the   decision  maker   by 
the  following  table  : 
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max  f.(x) 

max  f2(x) 

* 

f 

130.00 

250.00 

MAG (84%) 

108.00 

210.00 

£~ 

0.00 

0.00 

If  the  decision  maker  accepts  the  MAG  (84%),  it  would  be 
the  nondominated  preferred  solution.  Otherwise,  the 
decision  maker  can  set  his/her  own  initial  goals  according 
to  the  feasible  range  of  the  objective  function  values  or 
tradeoff  around  the  MAG.  Although  these  initial  goals  set 
according  to  these  reference  values  still  can  not  guarantee 
that  the  optimal  GP  solution  will  be  nondominated,  the 
results  should  be  much  closer  to  the  preferred  solution. 
The  decision  maker  also  can  save  a  lot  of  effort  in  deciding 
the  appropriate  initial  goals.  The  nondominance  of  the  GP 
solution  can  further  be  checked  by  the  techniques  Introduced 
in  the  next  section. 
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4.1.2  Techniques  for  Checking  the  Dominated  Solutions 

In  section  3.2.  we  conclude  that  the  GP  approach.  using 
preemptive  or  any  other  types  of  weighting  structure,  can 
lead  to  dominated,  suboptlmal  solutions  if  the  initial  goals 
are  too  conservative.  The  reason  is  that  the  GP  relies 
heavily  on  the  "principle  of  satisf icing"  and  the  decision 
makers'  perceptions  of  the  range  of  choice  and  feasibility. 
Consider  Figure  4.1.  which  presents  a  graphical  illustration 
of  this  shortcoming.  The  example  Is  in  objective  function 
space  (f  ,f  ).   The  region  where  feasible  solutions   lie   is 

ABCDEHIJA,  and  the  set  of  nondominated  points  lies  on  ABCDE. 
The  nondominated  set  is  unknown  to  the  decision  makers  when 
they   set   the   one-sided   goals   b    and   b    for   the   two 

objectives  E^      and   f 2.   respectively.     Since   (b  .b _)   i3 

feasible  in  Fq  in  Figure  4.1,  the  solution  to   the   GP   will 

give  total  deviations  of  zero;  i.e..  both  goals  are 
attained.  The  solution  is  inferior,  however.  since  it  is 
dominated  by  all  points  in  KLM.  Also,  point  C  is  an  extreme 
point  that  dominates  (b-.b.)   and   that   is   an  alternative 

optimum  of  the  GP  problem.  If  the  decision  makers  3tick 
with  their  original  goals,  then  they  will  settle  for  less 
than  they  should . 

Thus.  if  GP  is  chosen  as  the  desirable  solution 
technique,  this  inherent  disadvantage  relative  to  other 
approaches  must  be  eliminated.    That   is.   some   techniques 
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should  be  used  to  identify  the  dominance  of  the  GP  solutions 
and  also  generate  the  GP-nondominated  solutions  if  the 
original  GP  solution  is  dominated.  Techniques  for  this 
purpose  have  been  proposed  by  Hannan  [HI],  Cohon  [C6]  and 
Masud  and  Hwang  [Ml].  The  purpose  of  this  section  is  to 
introduce  these  techniques  and  identify  the  situations  to 
use  them. 


'1  "1 

Figure  4.1  Dominance  of  some  GP  solutions 
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The  full  expression  of  GP  model  after  normalization   can 
be  presented  by  : 

k 

Min  [    2     wrVi>  1  (4.8) 

r  =  l 
s.t. 

x  e  X 

f j(x)  +  (  <*]    "  <3j  )  hj  =  bj  ,      j  s  J 

^(S)  +  (  «j  -  A*   )  hx  «  ^  .     i  «  i 

fp(x)  +  (  a"  -  d+  )  hp  =  bp  .     p  e  p 
where 

f .(x)  :  the  benefit  goal  criteria,  j  s  J 

f,(X)  :  the  cost  goal  criteria,  i  m   I 

fp(£)  :  tne  two-sided  goal,  p  s  P 

d  ,  d   :  devlational  variables 

a^Cd)  :  achievement  function 

ar(d)  =  d~  ,       if  r  e  J 
a  (d)  =  d_  ,       if  r  e  I 


ar(d)  =  df  +  d*  ,  if  r  s   P 
b  :  initial  goals  set  by  decision  makers 
wf  :  relative  weights  among  objectives  and  can 

be  preemptive  or  nonpreemptive . 
h  :  the  coefficients  of  devlational  variables  of 
the  normalized  goal  functions. 
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hr  =  »crll  =  (cr-cr)1/2  =  Euclidean  norm 


The   following   methods  are   suggested   to   remedy   the 

shortcomings   of   dominated  GP   solutions.       Since    the 

normalization  is  considered,  the  original  formulations   will 
be  slightly  modified  here. 

Hannan'a  Procedures  : 

Hannan  [Hi]  extended  the  definition  of  a  GP  solution  to 
that  of  a  GP-nondominated  solution  and  demonstrated  how  to 
determine  when  dominated  GP  solution  occurs  and  how  to 
generate  the  set  of  GP-nondominated  solutions.   Suppose  that 

*'  =  (  < <    ><     *"      =      (   «][' <C   )•   d+*   ■   ( 

+  *       +* 

1  dn   }  is  an  optical  solution  to  (4.8).   Then,  define 

the  following  problem  to  check  the  dominance  of  GP   solution 


k 


Max 


s.t. 


I   <3i  +  v  .  u 


(4.9) 


r  =  l 
x  e  X 

',(«)  ♦  c  d"*  -  d+*  )  h.  .  b.  ♦  t .  .         3  e  j 
^Cx)  ♦   (  d"*  -  d+*  )  „x  =  bi  .  3i  ,         t  .  j 

V*>  +  <  %*  -  dP*  >  hP  ■  bP  • 
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t  >  0,  3  >  0, 

where  t  =  (  t% tk  )  and  s  =  (  3j «  ). 

We  are  now  ready  to  define  a  GP-nondominated  solution. 

Definition  of  GP-nondominated  Solution  : 

If  the  optimal  solution  to  (4.9)  is  zero,  the 
corresponding  optimal  solution  x  =  x*,  d"  =  d~*,  d+  =  d+* 
to  (4.8)  is  GP-nondominated. 

This  definition  implies  that  a  solution  to  (4.8)  is 
GP-nondominated  if  its  achievement  with  respect  to  two-sided 
goals  is  as  good  as  that  of  any  optimal  solution  to  (4.8), 
and  it  is  not  dominated  by  any  other  solutions  to  (4.8)  with 
respect  to  the  one-sided  goals. 

If  the  optimal  solution  to  (4.8)  is  dominated,  it  is 
necessary  to  generate  one  or  more  GP-nondominated  solutions 
for  consideration.  Hannan  proposed  the  following  MOLP 
problem  formulation  to  generate  the  GP-nondominated 
solutions  : 


«ax  f  .(x),     jsJ  (4  1Q) 

Min  fi(x),     i  «  i 
s.t. 

Ill 
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f1(x)   <  bx  +  hid+*    ,  j.  e  i 

V*>   '"   bp  "  V"p  +  hpdp*  ■      i  «  P 


where  x  =  x*.  d~  =  d~*.  d+  =  d+*  ia  the  dominated  optimal 
solution  to  (4.8).  (  Notice  that  the  GP  optimal  solution  13 
not  necessarily  a  nondomlnated  solution.  ) 

Therefore,  the  procedure  to  identify   the   dominated   GP 

solution  and  generate  the   GP-nondominated   solution   is   as 

follows.   First,  problem  (4.9)  should  be  solved  .   If  the  GP 

optimal  solution  is   identified   as   a   dominated   solution. 

(4.10)   should   be   solved.     The   nondomlnated    solutions 

generated  from  formulation  (4.10)  will  dominate  the  original 

GP  optimal  solution.    If  a  large  number  of  solutions  result 

from  (4.10).  this  is  an  indication  that  some  goals  have  been 

set  too  pessimistically.   This  can   be   remedied   by   either 

setting  more  demanding  goals  or  by  placing   weights   on   the 

objectives   in   (4.10).     These   weights   can    be    either 

preemptive  priorities  or  non-preemptive   weights.   but   they 

should   reflect   the   relative   value   of    extending    each 

objective  once  the  goals  have  been  satisfied  to  the   fullest 

extent. 

Masud  and  Hwang's  Procedure  : 

Another  method  proposed  by  Masud  and  Hwang  [Ml]  can  also 
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be  used  to  generate  the  nondominated  solutions.  This  method 
ensures  that  a  GP  solution  will  be  nondominated.  The 
formulation  with  slight  modification  is  as  follows  : 


n  [  (  I    wrV4>  )'  (  I    W2>  ]  1         (*•") 


Mi 
s.t. 


x  m   x 


f .j<*>  +  (  dj  ~  dj  >  hj  =  bj  '      3  e  J 

^(x)  +  <  dI  "  di  )  hx  =  b±  .      i  e  i 

fp(xj  +  (  d"  -  d+  )  hp  =  bp  .      p6p 

where  all  the  notations  are  the  same  as  (4.8)  and 

a«<d.)  =  <J   ,        if  s  «=  J 


,(d)  =  d3  '        if  s  e  I 


(d)  =  d"  +  d*  ,   if  s  e  P 


3 


This  procedure  is  very  similar  to  that  proposed  by 
Hannan.  The  first  preemptive  priority  in  this  formulation 
tries  to  satisfy  the  goals  first,  and  then  the  second 
preemptive  priority  guarantees  a  nondominated  solution  by 
optimizing  each  objective  with  weights  w  to  reflect  the 
relative  importance  of  extending  each  objective.  The 
difference  is  that  the  Masud  and  Hwang  procedure  produces 
a  single  solution  that  dominates   the   initial   GP   solution 
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whereas  the  Hannan  procedure  generates  all  solutions  that 
dominate  the  initial  GP  solution.  Masud  and  Hwang's 
procedure  is  convenient  and  easy  to  solve,  while  Hannan's 
procedure  of  using  MOLP  fomulatlon  is  more  complicated  and 
it3  computer  package  is  not  aa  commom  as  LP. 

Below,  we  use  the  same  numerical  example  which  has  been 
used  in  section  3.2  to  illustrate  the  procedures  of  checking 
the  dominated  solutions  and  generating  the  nondominated 
solutions. 


Numerical  Example  4.2  :  (  Hardee  Scheduling  problem  ) 

The  Hardee  Company  makes  two  kinds  of  dolls.  Doll  A  is 
a  high  quality  one  and  Doll  B  is  of  lower  quality.  The 
respective  profits  are  30.4  and  30.3  per  doll.  Each  Doll  A 
requires  twice  as  much  time  as  a  Doll  B,  and  if  all  dolls 
were  of  type  B,  the  company  could  make  500  per  day.  The 
supply  of  material  is  sufficient  for  only  400  dolls  per  day 
(  both  A  and  B  combined  ).  The  problem  assumes  that  all  the 
dolls  for  type  A  and  type  B  the  factory  can  make  could  be 
sold,  and  that  the  best  customer  of  the  company  wishes  to 
have  as  many  as  possible  of  type  A  doll.  The  manager 
realizes  that  two  objectives  :  (1)  the  maximization  of 
profit  and  (2)  the  maximum  production  of  Doll  A.  should  be 
considered  in  the  scheduling  the  production. 

By  using  the  GP  approach  without   first   exploring   the 
feasible  solution  space,  the  manager  establishes  the   goals 
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for  the  two  objectives  as  follows  : 

Gl  :  achieve  the  profit  of  at  least  390. 
G2  :  produce  at  least  180  units  of  Doll  A. 

It  is  assumed  the  relative  importance  of  the  two 
objectives  are  the  same,  and  vector  normalization  is 
applied.   Then,  the  GP  formulation  (4.8)  for  this  problem  is 


min 


[  «l   +  *a  ] 


0.4x1  +  0.3x2  +  h^d"  -  d*)  =  90 
Xl  +  Vd2  "  d2>  =  180 


x1  +  x2  <  400 


2xx  +  x2  <  500 


X,  d~.  d+  >  0 


where    h1    =    ||cr||    =    /o.42+0.32         =0.5      and   h      >    1. 

The  solution  of  this  GP  model  is  : 

x   =  (  225,  0  ),   d~*  =  (  0,  0  ),   d+*  =  (  0.  45  )  and 
I  =  (  90.  225  ). 

This  solution  is  point  F  in  Figures  3.2  and  3.3. 


Now,  we  are  ready  to  check   the   dominance   of   this   GP 
solution. 
Hannan'3  procedure  : 

Step  1^  :  Check  the  dominance  of  GP  optimal  solution. 
[  by  using  the  formulation  (4.9)  ] 


max  U  =  (  tx  +  t2  ) 


s.t. 


0.4x1  +  0.3x2  +  0.5  (  0  -  0  )  =  90  +  t 
x1+l(0-45)=  180  +  t 

The  solution  to  this  problem  is  : 

x  =  (  250,  0  ),  t  =  (  10,  25  ),  and   U  =  35. 

Because  U   does   not   equal   to   zero,   the   GP   optimal 

* 
solution  x   is  dominated  solution.   So,  the  next  step  is   to 

generate  the  GP-nondominated  solutions. 


Step  2  :  Generate  GP-nondominated  solutions. 
[  by  using  the  formulation  (4.10)  ] 


max  0.4x   +  0.3x 


s.t. 

x  e   X 

O^Xj^  +  0.3x   >  90  -  0.5  (  0  ) 

x1  >  180  -1(0) 
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The  GP-nondominated  solutions  to  this  problem  are  : 

CD  X  =  <  225,  50  ).  f  =  (  105,  225  ).  i.e..  point  H 
in  Figures  3.2  and  3.3. 

(2)  x  =  (  250,  0  ).  f  =  (  100,  250  ),  i.e..  point  E 
in  Figures  3.2  and  3.3. 

Both  of  the  nondominated  solutions  dominate  the  original   GP 
solution  which  is  x  =  (  225.  0  )  and  f  =  (  90,  225  ). 


Masud  and  Hwang's  procedure  :  [  formulation  (4.11)  ] 
Ln   [  (  d^+d^  ),  (  -d+-d*  )  J 


mir 


i.t. 

x  m  X 


0.4Xl  +  0.3x2  t  hj  (  dj  -  d|  )  =  90 
x1  +  h2  (  d~  -  d*  )  =  180 

where  hj  =  ||cC||  =  /o.42+0.32   =  0.5  and  h£  =  1. 
The  nondominated  solution  generated  by  the  this  procedure  is 


x  =  (  250.  0  ).  f  =  (  loo.  250  ).   i.e..   point   E   in 
Figures  3.2  and  3.3.. 
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In  general,  a  set  of  goals  will  lead  to  a  nondominated 
GP  solution  only  if  the  goals  are  in  the  nondominated  set  or 
beyond  the  feasible  region.  Thus,  whenever  a  GP  solution 
yields  zero  deviations,  the  analyst  should  auspect  that  the 
set  of  goals  has  led  to  dominated  solution.  The  methods 
described  in  this  section  are  Important  because  it  provides 
an  avenue  for  dealing  with  dominated  GP  solutions.  It  would 
be  the  analysts'  responsibility  to  check  the  nondominance  of 
the  GP  solution  and  generate  the  possible  GP-nondominated 
solutions  for  consideration.  The  techniques  introduced  in 
this  section  are  suggested. 
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4.1.3  Use  of  Nonpreemptive  Weighting  Approach 

The  preemptive  weighting  scheme  allows  the  decision 
maker  to  prioritize  goals  on  an  "ordinal"  basis.  From  a 
practical  point  of  view,  it  seems  relatively  easy  for  the 
decision  maker  to  rank  the  importance  of  multiple  goals. 
However,  the  improper  use  of  this  weighting  scheme  may  cause 
more  harm  than  good.  In  section  3.3  and  3.4,  we  discussed 
the  shortcomings  of  preemptive  GP.  The  weakness  of 
preemptive  GP  emerges  in  three  ways  :  (1)  The  preemptive 
approach  is  not  compatible  with  the  nonpreemptive  approach 
and  the  decision  maker  runs  the  risk  of  not  obtaining  the 
preferred  solution.  (11)  The  preemptive  models  may  totally 
ignore  lower  priority  objectives.  (ill)  The  preemptive 
models  disallow  tradeoffs  of  small  losses  in  high  priority 
objectives  for  large  gains  in  low  priority  objectives. 
Thus,  the  preemptive  GP  is  not  recommended  in  solving  the  GP 
problems  except  that  the  goals  really  do  have  preemptive 
importance  in  relationship  to  one  another. 

If  GP  has  to  be  used,  nonpreemptive  GP  would  be 
better  than  preemptive  GP.  GP  with  nonpreemptive  weights 
Is  merely  a  LP  model  and  can  be  solved  easily  and 
efficiently.  The  extensions  such  as  duality,  sensitivity 
analysis,  and  parametric  programming  can  be  handled  in  the 
same  manner  as  in  LP  problems.  The  only  difficulty  of  using 
nonpreemptive  GP  is  the  determination  of  a  set  of  accurate 
weights  which  reflect  the  priorities  of  the  various  goals. 
The  purpose  of  this  section   is   to   suggest   the   auxiliary 
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techniques  to  help  the   decision   makers   in   assessing   the 
weights. 

In  practice,  the  selection  of  appropriate  weights  does 
not  appear  to  cause  GP  problem  formulators  too  much 
difficulty  as  it  does  to  the  formulators  in  using  other 
utility  preference  methods.  The  presumption  of  this  is  that 
if  GP  is  chosen  as  the  desirable  solution  technique,  the 
problem  usually  tend  to  have  explicit  (  or  implicit  ) 
reasons  for  the  scale,  weights  and  trade-offs  between  goals, 
and  the  decision  makers  should  be  quite  familiar  with  the 
objectives  so  that  he/she  can  set  the  appropriate  initial 
goals  and  is  able  to  provide  trade-off  information  with 
which  an  acceptable  set  of  weights  can  be  determined, 
otherwise  GP  should  not  be  used.  However,  as  the  goals 
often  represent  different  units  and  the  decision  maker's 
judgements  can  not  be  so  accurate  and  consistent,  the  direct 
weighting  of  the  GP  objectives  is  questionable.  Therefore, 
the  normalizing  and  weighting  procedure  must  be  applied  to 
reduce  the  burden  of  decision  maker.  The  eigenvector  method 
is  suggested  for  weighting  the  relative  importance  of  goals 
in  this  section. 

Normalization  : 

Normalization  provides  the  comparable  scale  for  the 
objectives  and  presents  a  convenient  way  to  correctly  set 
the  desired  weights.  Thus,  the  goal  constraint  function 
f(x)  should  be  normalized  before   determining   the   relative 
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weights  to  assign  to  d  and  d  in  the  objective  function. 
The  objective  function  relative  weights  would  be  determined 
on  the  basis  of  how  important  it  is  to  achieve  each  goal 
relative  to  the  Importance  of  achieving  the  other  goals. 
This  is  an  important  consideration  in  developing  models 
which  are  as  representative  as  possible  of  the  concerns  of 
the  decision  maker.  Two  different  schemes  of  scale 
normalization  have  been  3tated  in  section  2.3.4. 

Eigenvector  method  : 

Another  auxiliary  technique  suggested  to  solve  the 
weighting  problem  in  nonpreemptlve  GP  is  the  eigenvector 
method.  This  process  is  intimately  connected  with  the  idea 
of  consistency  of  thought.  Its  mathematical  foundations  are 
simple.  Satty  [SI]  introduced  this  method  by  using 
the  eigenvector  principle  in  a  positive  pairwlse  comparison 
matrix.  By  using  the  eigenvector  method.  decision  maker 
just  needs  to  compare  the  relative  importance  of  two 
objectives  (  or  goals  )  each  time  and  the  pairwlse 
comparison  matrix  will  be  constructed  according  to  decision 
maker's  judgements.  If  the  number  of  objectives  is  n,  the 
number  of  judgements  is  rC2  =  n(n-l)/2.  The  principle  of 
eigenvector  is  stated  below.  Let  the  pairwlse  comparison 
matrix  A  be 
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A   = 


11 
»21 


12 
*22 


In 

'2a 


nl 


Thl3   is   a   "reciprocal   matrix"   which   has   all   positive 
elements  and  has  reciprocal  property: 

=  1  / 


U 


ji 


(4.12) 


and  consistency  property  : 


lj 


ik 


jk 


(4.13) 


Multiplying  A  by  w  =  (w  ,  w  ,, 


wn)   yields 
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"wl ' 

■    n 

"Wl    ' 

W2 

W2 

w 
n 

w 

n 

or    (  A  -  nl  )w 


(4.14) 


Due  to  the  consistency  property  of  eq.  (4.13).  the  system  of 
homogeneous  linear  equations,  eq.  (4.14),   has   only   trlval 

solutions,  that  Is,  all  the  eigenvalues  X.  .  i  =  1   n, 

of  matrix  A  are  zero  except  one  and  the  value   of   the   only 
nonzero  eigenvalue  i3  equal  to  n. 

In  practice,  the  precise  values  of  w./w.  are  unknown  and 
must  be  estimated  by  decision  maker.  In  other  words,  human 
judgements  can  not  be  so  accurate  that  eq.  (4.13)  be 
satisfied  completely.  We  know  that  in  a  positive  reciprocal 
matrix,  small  perturbations  in  the  coefficients  imply  small 
perturbations  in  the  eigenvalues.  Hence  the  eigenvector  is 
insensitive  to  small  changes.  If  we  define  A'  as  the 
decision  maker's  estimate  of  A  and  w"  i3  corresponding  to 
A' ,  then 


A'w'  =  X   w' 
—     max— 


(4.15) 


where  Xmax  is  the  largest  eigenvalue  of  A'.    The   weight! 


ng 


ioa 


vector  w'  can  be  obtained  by  solving  the  system  of  linear 
equations,  eq.  (4.15). 

For  assessing  the  scale  ratio  w  /w  . ,  Saaty  [SI]  gives  an 

intensity  scale  of  Importance  for  activities  and  has  broken 
down  the  importance  ranks  as  shown  in  the  scale  in  Table 
4.1. 

Table  4.1   Intensity  of  Importance  Scale 


Intensity  of 
Importance 


Definition 


Explanation 


,4.6.8 


Equal  importance 


Weak  importance  of 
one  over  another 


Essential  or  strong 
importance 


Demonstrated 
importance 


Absolute  importance 


Intermediate  values 
between  the  two  ad- 
jacent judgements 


Two  activities  contribute 
equally  to  the  objective 

Experience  and  judgement 
slightly  favor  one 
activity  over  another 

Experience  and  judgement 
strongly  favor  one 
activity  over  another 

An  activity  is  strongly 
favored  and  its  dominance 
is  demonstrated  in  practice 

The  evidence  favoring  one 
activity  over  another  Is 
of  the  highest  possible 
order  of  affirmation 

When  compromise  is  needed 


This  weighting  process  is  illustrated  by   the   following 
nutrition  problem  example  [C4]. 
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Numerical  Example  4.3 :  (  Nutrition  Problem  ) 

The  nutrition  problem  is  to  determine  the  quantities  of 
six  foods  that  should  be  eaten  to  balance  nutritional 
requirements.   The  three  objectives  of  this  problem  are  : 

Goal  1  :  achieve  carbohydrate  intake  at  least  450  cal. 
Goal  2  •.  achieve  cholesterol  intake  at  most  28  units. 
Goal  3  :  satisfy  the  cost  goal  not  exceeding  9  3.0. 

In  this  example,  we  put  the  stress  on  the  weighting  of 
these  three  goals  by  eigenvector  method  and  pay  no  attention 
to  the  problem  formulation  and  solution  process. 

Suppose  the  positive  pairwise  comparison  matrix  is 
estimated  by  decision  maker  as  : 


A  = 


Gl 

G2 

G3 

1 

1/3 

1/2   ' 

Gl 

3 

1 

3 

G2 

2 

1/3 

1 

G3 

then  set  the  determinant  of  (A-\I)  as  zero.   That  is 


det(A-XI)  = 


-A 

1/3 

1/2 

3 

1-X 

3 

2 

1/3 

1-X 

The  largest  eigenvalue  of  A,  X    .  is  3.0536 

max 


and  we  have 
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-2.0536       1/3       1/2 
3         -2.0536       3 
2  1/3      -2.0536 


=  0 


The  solution  of  the  homogeneous  system  of   linear   equations 

..T 


gives  (  recall  that  w  +  w  +  w 


1  )  ■•  w'  =  (  0.1571, 
0.5936.  0.2493  ).  This  Is  the  weighting  vector  for  the 
three  goals. 


Consistency  in  judging  the  relative  importance  of  two 
goals  is  necessary  to  get  valid  results.  Satty  [SI] 
suggested  a  measure  of  the  consistency  or  reliability  of 
judgements  as 


C.I. 


X    -  n 
max 


and  the  average   consistencies   for   different-order   random 
matrices  are  : 


Size  of  matrix 


Random  Consistency 


1 
2 

3 

4 

5 
6 

7 

a 

9 
10 


.00 

.00 

.58 

.90 

1.12 

1.24 

1.32 

1.41 

1.45 


If  we  divide  C.I.  by  the  random  consistency  number  for  the 
same  size  matrix,  we  obtain  the  consistency  ratio  (  C.R.). 
The  value  of  C.R.  should  be  around  10  percent  or  less  to  be 
acceptable.  If  it  is  not  acceptable,  we  may  need  to  revise 
our  judgements.   In  our  example 


„  T       3.0536  -  3       „  „„ 
C.I.  =   - ; =  0.0268 


„  „         0.0268 

C.R.  =   g-gg =4.6%    (  acceptable  ) 


This  suggested  weighting  procedure  can  really  help  the 
decision  maker  make  a  consistent  decision  about  the 
importance  of  the  goals  (  or  objectives  ).  But  consistency 
does  not  mean  getting  an  answer  closer  to  the  preferred 
solution  of  decision  maker.  Therefore,  when  the 
nonpreeraptlve  GP  is  used,  the  results  must  Include  a 
sensitivity  analysis  and  multiple  runs  with  different  sets 
of  weights,  so  that  a  reasonable  number  of  compromise 
solutions  are  available  to  the  decision  maker  and  the 
decision  maker  can  have  more  confidence  in  the  selected 
preferred  solution. 
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4.1.4  Interactive  Approach  of  GP  —  IS6P  II 

Interactive  techniques  for  solving  multiple  objective 
problems  are  based  on  a  search  process  which  consists  of 
periodic  feedback  from  the  decision  maker,  which  serves  to 
guide  the  direction  of  the  search.  Each  cycle  in  the 
process  contains  an  optimization  phase  and  a  decision  maker 
evaluation  phase.  The  true  preference  function  of  the 
decision  maker  remains  unspecified,  but  the  direction  of  the 
decision  maker's  preference  is  identified  at  each  stage  by 
presenting  the  decision  maker  with  a  trial  solution.  The 
advantages  of  an  -interactive"  approach  are  that  they 
provide  for  the  decision  maker  a  learning  process  about  the 
problem  and  make  allowance  for  psychological  convergence  for 
the  decision  maker.  Interactive  GP  techniques  are  defined 
as  those  interactive  techniques  in  which  each  cycle  contains 
the  development  of  a  set  of  new  goals  specified  by  the 
decision  maker. 

In  this  chapter.  some  remedial  methods  have  been 
suggested  to  make  the  shortcomings  of  standard  GP  less 
severe,  but  none  of  them  are  as  effective  as  the  interactive 
GP  approach.  GP  with  interactive  approaches  does  not  have 
many  of  the  practical  and  theoretical  shortcomings  of  the 
standard  GP.  The  fact  that  there  is  a  formalized  procedure 
for  generating  new  solutions  in  the  interactive  GP  methods 
mitigates  many  of  the  criticisms  of  GP.  In  many  interactive 
GP  approaches,  preemptive  GP  is  not  used  and  the  setting  of 
the  initial  goals  is   not   as   objectionable   because   these 
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goals  are  inevitably  altered  during  the  interactive  process 
as  the  decision  maker  learns  more  about  the  decision  space. 
Also,  the  non-dominated  GP  solution  is  guaranteed. 

The  Interactive  approach  of  GP  is  Intuitively  appealing 
because  the  features  of  GP  can  enhance  the  efficiency  of 
finding  the  best  compromise  solution, and  the  decision  maker 
is  more  cognitively  oriented  toward  achievement  for  the 
objectives  than  toward  Indirect  measures  such  as  weights. 
Several  MODM  methods  have  been  created  in  attempts  to 
provide  a  link  between  GP  and  interactive  approaches  such  as 
:  Interactive  Goal  Programming  (IGP,  Dyer.  1972)  [D2],  Goal 
Programming  STEM  (GPSTEM,  Fichefet,  1976)  [Fl],  Sequential 
Information  Generator  for  Multiple  Objective  Problems 
(SIGMOP,  Monarchi.  et  al.,  1975)  [M2],  Interactive  Multiple 
Goal  Programming  (IMGP,  Spronk,  1981)  [S5],  Interactive 
Sequential  Goal  Programming.  (ISGP,  Masud  and  Hwang,  1981) 
[Ml],  and  Interactive  Sequential  Goal  Programming  II  (ISGP 
II,  Hwang  and  Masud,  1987)  [H7]. 

In  this  section,  the  ISGP  II  is  recommended 
to  solve  the  GP  problems  because  it  combines  the  attractive 
features  of  Goal  Programming  and  other  interactive  MODM 
methods.  ISGP  II  is  based  on  the  Implicit  assumption  that 
the  decision  maker  can  adjust  his/her  desired  goals  through 
an  iterative  learning  process  based  on  information  in  a  set 
of  solutions.  The  information  provided  in  each  iteration 
contains  a  current  best  compromise  solution,  called 
Principal  Solution  (PS),  and  a  set   of   Auxiliary   Solutions 
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(AS)  which  are  compromise  solutions  obtained  if  each  of  the 
goals  set  by  decision  maker  are  satisfied  serially.  The 
Positive  Ideal  Solution  (PIS),  Negative  Ideal  Solution 
(NIS),  Maxlmuim  Achievable  Rate  (MAR),  and  Maximum 
Achievable  Goal  (MAG)  for  the  multiple  objectives  are  also 
used  as  the  reference  values  to  help  the  decision  maker  set 
the  goal  levels.  The  PIS.  NIS,  MAG,  and  MAR  are  determined 
prior  to  having  the  decision  maker  set  the  initial  goals. 

ISGP  II  explores,  selectively,  the  nondomlnated 
solutions  of  MODM  problem  in  an  interactive  way.  Any 
iteration  r  consists  of  two  phases  :  calculation  phase  and 
evaluation  phase.  The  PS  and  AS  are  obtained  in  the 
calculation  phase  and  the  evaluation  phase  consists  of  the 
decision  maker's  indication  of  his/her  preference  judgement 
about  these  solutions  in  the  form  of  new  desirable  goal 
levels.  With  this  new  Information,  the  process  goes  back  to 
the  calculation  phase  of  the  (r+l)th  iteration.  If  in  any 
iteration  the  decision  maker  finds  that  either  the  PS  or  any 
of  the  AS  is  a  satisfactory  solution,  the  algorithm 
terminates.  If  the  decision  maker  provides  his/her 
preference  information  in  a  consistent  manner,  then  the 
number  of  iterations  necessary  is  expected  to  be  small. 

The  algorithms  of  ISGP  II  are  described  below  : 
Problem  Formulation  : 


The  MODM  problem  is  : 

Max/Min   [  f^x),  f2(x) fk(x)  ] 

s.t.      x  e  X 

where   £\(x)  :  the  benefit  objective,  j  s  J 
fi(X)  :  the  co3t  objective,  is  I 

Algorithm  : 

(  The  algorithms  of  PIS,  NIS,   MAR,   and   MAG   have   been 
stated  in  section  4.1.1  and  are  not  repeated  here  again.  ) 

Step.  0  :  Determine  the  Positive   Ideal   Solution   (PIS),   f*, 
and  Negative  Ideal  Solution  (NIS),  f"  : 
(  see  section  4.1.1  ) 

Step  1  :  Set  initial  goals. 

Step  JUJ.  :  Obtain  a  set  of  Maximum  Achievable  Goals  (MAG) 
and  a  Maximum  Achievable  Rate  (MAR),  B%.  (  The  MAG  is  a 
feasible  solution  which  has  the  same  rate.  MAR,  for  all 
objectives.  ) 

(  also  stated  in  section  4.1.1  ) 


Step  1.2  :  Display  ISGP  II  Inlti 


al  Table. 
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max/mm    max/min  max/min 


"kl 


MAG(B%)      b*        b*  .  .  b* 

12  dk 

<-  'I         rt  ...  f- 


The  initial  table  is  presented  to  the  DM.  If  the  DM  is 
satisfied  with  the  MAG  then  take  it  as  the  preferred 
solution  and  go  to  step  8.  Otherwise  go  to  step  1.3. 

steP  1-3  :  Set  initial  goals. 

The  DM  set  his/her  initial  goals  at  b  referring  to  the 
initial  table.   The  initial  goals,  b.  should  satisfy 


f ~i      <biS      f i    ■       1    *J 


Steja  2  :  Solve  for  PS  (Principal  Solution). 

The  problem  is  to  find  a  set  of  x  so  as  to 


(4.16) 


mln  a 

ps 


r  =  l 
a-t.    x  e  X 


[  I     Wr  dr   ]  (4.17) 


*j(5>  +  (dj  "  dj>  «  f  ■  "  V  -  bj 


J  s  J 
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fi(X)   -   (di   -  dj)   (   f"   -  f*)   *  ^   ,   i  «  J 

0  -  dr  '  dr  "  1 

where  «r  is  the  relative  importance   of   each   objective 
33igned  by  the  DM. 


Assume  the  solution  of  (4.17)  is  :  x°,  f°   a 

-    ps 


Ster.  2.1  :  Check  if  f   is  a  set  of  nondoralnated  solutions  : 

If  apg  =  °»  then  f   is  a  set  of  dominated  solutions,  and 
go  to  step  2.2. 

If  aps  "   °'  tnen  take  1°  a3  the  principal  solution  and  go 
to  3tep  3. 

Step  2.2  :  Improvement  of  f°. 

We  increase  b^  and  decrease  bx  by  z%.  recalculate 
(4.17),  and  check  the  value  of  apg.  The  procedure  repeats 
until  an  apg  *  0  is  found.   z%  is  recommended  to  be   set   at 

5%.      The  goal  setting  rules  in  (4.16)  must  be  satisfied  when 
modifying  b .  and  b  . 


Ster.  3  :  Solve  for  k  sets  of  ASs  (Auxiliary  Solution). 

Each  set  of  auxiliary  solutions  is  obtained  by  setting 
one  individual  goal  as  a  constraint  and  by  minimizing  the 
summarized  weighted  deviations  from  other  goals. 
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The  problem  is  to  find  a  3et  of  x  30  as  to 

k 


min  au   =  [  I     wr  d" 


(4.18) 


r  =  l 
r*u 


3.t.    x  e  X 

fjCX)  +  (d"  -  d+)  (  f*  -  f",  =  b_.  ,(for  all  j,  3  „  u) 

fi<*)  "  (d^  -  <**)  (  f~  -  f*)  =  b±  .(for  all  i.  i  *  u) 

fuC2)  i  bu   (  If  u  i  J  ) 

fu(-)  s   bu   (  if  u  s  I  ) 


0  <  <T  .  ac  <  1 

Assume  the  solutions  of  (4.18)  are 


C  ,  £  .  a  ,  U  a  1,2 it 


Stee  3_;_l  :  Check  if  fu  i3  a  set  of  nondominated  solutions. 
If  a   =0,  then  fu  is  a  set  of   dominated   solutions   and 
improve  the  b.    and  b±  as  the  same  way  as   step 

2.2  until  an  a  *   0  is  found. 
If  au  x   0,  then  take  fu  as  the  uth  auxiliary  solution   and 
go  back  to   step  3   for   the  (u+l)th   auxiliary 
solution  until  all  k  gets  of  auxiliary  solutions 
are  obtained,  u  =  1,2 k.   Go  to  step  4. 


119 


Step  4  :  Form  ISGP  II-Tradeoff  Table. 

The  tradeoff  table  ia  given  in  Table  4.2  to  show  the   DM 
all  the  nondominated  solutions  PS,  and  ASs. 


Table  4.2   ISGP  II   Tradeoff  Table 


max/rain    max/min                  max/min 

eAcx)         f2(x)           •  •   •                fk(x) 

*         fi      f2       •■•  •        *;; 

MAG(B%)          b*        b*          ...           b* 

i"                    fi            f2                •   •  •                 '» 
----- !i i_        •  •  • 

Iteration  # 

g°al            bl        b2          •  •  •           \ 

^        f°          ...           f° 
f2          •  •  "           *jj 

fi                               •   •  •                  e* 

£i             f* 

Step  ji  :  Evaluation  phase. 


Step.  5^  :  Ask  the  DM  :  -  Is  f°  or  any  fu  satisfactory  ? 

If  yes.  then  go  to  step  8,  otherwise  go  to  step  5.2. 
Step.  5^2  :  Ask  the  DM  to  specify  new  goals.  b,  from  the 
information  in  the  tradeoff  table. 


Step  6  :  Calculate  the  deviations  of  the  new  goals  from  PS. 


r.  -  f°  -  b 


i    .    for  all  j  =  {  j  |  b   <  f°  } 


f j  =  bj  "  f j  '  f°C  a11  j  =  {  j  '  bj  >    f j  } 


"i 


b±  -  £L    .    for  all  i  =  {  i  |  b±  <  f°  } 


fi  "  fi  "  bl  '  for  all  1  =  {  i  |  bx  <  f°  } 

Step  2  !  Consistency  check  for  the  specified  b. 

For  all  r  =  1.2 ,  k.  check  to  see  if  there  is  at  least 

one  f   >  0  when  one  or  more  f   >  0. 

If  yes,  then  go  back  to  Step  2  and  start  the  next 
iteration,  otherwise  go  to  Step  5.2.  Since  f°  is  a 
nondominated  solution,  this  consistency  check  is  necessary 
to  make  sure  that  if  the  DM  wants  an  improvement  in   one   of 

the  objective  values  from  f°,  s/he   is   ready   to   accept   a 
detriment  in  at  least  one  other  objective  value. 

Step  8.   Stop. 

Print  the  preferred  solution  x  and  f   and   terminate   the 
algorithm. 

Numerical  Example  4.4   (  continued  from  Example  4.1  ) 

The  same  numerical  example  of  section  4.1.1  is  used  to 
illustrate  the  solution  procedures  of  ISGP  II.  Since  the 
PIS,  NIS.  MAG,  and  MAR  have  been  solved  in  Example   4.1,   we 
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start  the  procedure  from  Step  1.2  Display   ISGP   II   Initial 
Table. 

Step  1.2  :  Display  ISGP  II  Initial  Table. 


max/mln 

* 

f 

130.00 

MAG(84%) 

108.00 

f 

0.00 

w 

1.00 

max /i 

nin 

e2(i 

£) 

250 

.00 

210 

.00 

0 

.00 

1 

.00 

We  assume  the  DM  does  not  accept   the   MAG (84%)   as   the 
preferred  solution. 

Step  1.3  :  Set  initial  goals  b. 

Assume  the  DM  sets  the  initial  goals  b  at  90%  of  |f* -  f " | 

,  i.e.,  b  =  (117.00,225.00).  referring  to  the  ISGP  II  initial 
table. 

Iteration  No.  1^  : 

Step  1  :  Solve  for  PS  (Principal  Solution). 

min   apg    =    {    d"    +   d^    >  (4.19, 

s.t.         x    e  X 


f1(x)    +    130    (d~    -    d*)    =    i>1    =    117 


f2(x)    +    250    (d~    -    d+)     =    b2    =    225 


0   <  ax    .    d2   <   1 


The  solution  ia  given  as  : 

x  =  (225.00,  50.00) 

£°  =  C 105. 00,  225.00) 

a    *   0 
pa 


Step   2.1  : 

As  apa  *   0  implies  that  the  solution  f°  is  a  nondominated 
3o!ution  and  la  accepted  as  the  PS.   Go  to  Step  3. 


Step  3  s  Solve  for  ASs  (  Auxiliary  Solutions  ). 

The  first  auxiliary  aolution  ia  obtained  by  solving 

mln  a   =  {  d2  }  (4.20) 

a.t.   x  m   X 

f1(x)  >  b1  =  117 

f2(x)  +  250  (d~  -  d*)  =  b2  =  225 

0  <  d"  ,d+  <  1 


The  aolution  ia  given  as  : 

x  =  (165.00,  170.00) 


f  =    (117.00,  165.00) 


.x  m   o 
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Step  3.1  : 

The  answer  a  *   0  implies  f1  is  a  nondominated   solution 
and  is  accepted  as  the  first  AS.   Go  back  to  Step  3  for   the 


second  AS  set. 

Steg  3  :  Solve  for  ASs  (  Auxiliary  Solutions  ). 

The  second  auxiliary  solution  is  obtained  by  solving 
min  a2  =  {  d"  }  (421) 

a.t.   x  e  X 


f1(x)  +  130  (d~  -  d+)  =  bx  =  117 


f2(x)  >  b2  =  225 

0  <  d~  .d*   <  1 


The  solution  is  given  as  s 
x  =  (225.00,  50.00) 

2 
f   =  (105.00,  225.00) 

a   *  0 


Step  3.1  : 

Since  a  *   0  implies  f2  is  a  nondominated   solution   and 
is  accepted  as  the  second  AS.   Go  to  Step  4. 

Step.  4  :  Form  ISGP  II  Tradeoff  Table. 
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max/min  max/min 

^(x)  f2Cx) 


f.  130.00  250.00 

MAG(84%)  108.00  210.00 

f  0.00  0.00 

H  1-00  1.00 

Iteration  1 

Goal  117.00  225.00 

PS  105.00  225.00 

AS  1  117.00  165.00 

AS  2  105.00  225.00 


This  tradeoff  table  le  presented  to  DM  for  evaluation. 

Step  5    :    Evaluation  phase. 

We  assume  the  DM  is  satisfied  with  none  of  the 
solutions.  We  proceed  to  Step  5.2  for  setting  another  set 
of  goals  b. 

Step  5.2  :  Set  a  set  of  new  goals  b. 

Since  none  of  the  PS  and  ASg  are  satisfactory  to  the  DM, 
a  set  of  new  goals  b  is  set  at  (110.00.215.00)  by  the  DM. 
Steg  6  :  Calculate  the  deviations  of  the  new  goals  from  PS. 

The  new  goals,  b  =  (110.00.  215.00)  are  compared  with  f° 
=  (105.00,  225.00)  of  the  first  Iteration. 

f1  =  110.00  -  105.00  =  5.00   (  b   >  f°  ) 

f   =  225.00  -  215.00  =  10.00  (  b   <  f°  1 

*  2     2  ' 
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Step  7  :  Consistency  check  on  new  b. 

Since  there  are  one  f"  >  o  and  f  +  >  0,  the  new  goals  b 
is  accepted  as  the  initial  goals  for  the  second  iteration. 
So.  start  from  Step  2  again. 


Iteration  2   : 

Step  2  :  Solve  for  PS. 

The  problem  of  (4.19)  is  solved  again  with  b    =   HO. 00 
and  t>2    =  215.00.   The  solution  is  given  as  : 

x  =  (215.00.  70.00) 
f°  =  (107.00.  215.00) 
apa  *   0 

Ste£  U  :  The  answer  aps  x  0  implies  that  f°  is  a  set  of 
nondominated  solutions  and  is  accepted  as  the  principal 
solution.   Go  to  Step  3  for  ASs. 

Step  3  :  Solve  for  ASs. 

Similar  to  that  of  iteration  1.   We  solve  (4.20)  with  b 
=  110.00  and  b2  =  215.00.   The  solution  is  given  as  : 

X  =  (200.00.  100.00) 


f1  =  (110.00,  200.00) 

a  "   0   (  ■•-  f   Is  accepted  as  the  first  AS  ) 


For  the  second  AS.  we  solve  (4.21)  with  b   =  110.00   and 


»2  =  215.00.   The  solution  is  given  as 


x  =  (215.00,  70.00) 

f2  =  (107.00.  215.00) 

2  2 

a  *   0   (  .-.  f   is  accepted  as  the  second  AS  ) 

Go  to  Step  4 . 


Ster,  4  :  Form  the  ISGP  II  Tradeoff  Table  for  iteration  2. 

max/rain  max/nln 
fjCx)  f2(x) 


I  130.00  250.00 

MAG(84%)  108.00  210.00 

I  0.00  0.00 

w  1.00  1.00 

Iteration  2 

Goal  110.00  215.00 

ps  107.00  215.00 

As  1  110.00  200.00 

AS  2  107.00  215.00 


This   tradeoff   table   is   presented   to   the    DM    for 
evaluation. 
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Step  5  :  Evaluation  phase. 

We  assume  the  DM  is  satisfied  with  the  first  auxiliary 
solution  (AS  1)  and  selects  it  as  the  preferred  solution. 
Go  to  Step  a. 

Step  8  :  Stop. 

The  final  solution  is  printed  as 

x  =  (200. 00,  100.00) 

f  =  (110.00.  200.00) 
and  the  algorithm  terminates  here. 

Goal  Programming  incorporate  the  interactive  approach 
can  remedy  most  of  the  shortcomings  of  the  standard  GP.  In 
particular,  most  interactive  GP  approaches  do  not  involve 
either  the  use  of  preemptive  priorities  or  of  a  priori 
Archimedean  weights.  This  relieves  the  decision  maker  of  a 
substantial  burden.  The  only  a  priori  requirement  of  the 
decision  maker  for  most  of  these  approaches  is  the 
specification  of  a  set  of  trial  goals.  Also,  in  these 
approaches  the  decision  maker  is  generally  provided  in 
advance  of  the  best  and  worst  possible  attainments  of  the 
objectives  as  an  aid  in  specifying  these  goals.  The  use  of 
these  goals  in  each  iteration  guides  the  direction  toward 
the  decision  maker's  best  compromise  solution.  Therefore, 
if  the  decision  maker  is  willing  to  invest  time  in  the 
interactive  process,  he/3he  can  always  find  the  feasible 
best  compromise  solution. 
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4.2  Identification  of  Using  GP  and  Alternatives 

To  solve  the  MODM  problems,  a  lot  of  methods  have  been 
created.  Each  method  has  different  assumptions  and 
different  Information  requirements.  In  Chapter  2,  Figure 
2.3,  many  of  the  MODM  methods  were  listed.  In  Chapter  3, 
many  of  the  shortcomings  were  also  discussed.  Nevertheless, 
GP  has  been  almost  the  most  popular  MODM  method  today.  The 
reasons  might  be  partially  due  to  its  simplicity.  easy 
usage,  and  its  shortcomings  not  being  well-known.  In 
addition,  a  lack  of  guidelines  in  helping  users  to  select 
the  appropriate  MODM  algorithm  might  be  another  reason. 
There  must  exist  some  situations  in  which  GP  should  not  have 
been  used.  Therefore,  there  is  a  need  to  investigate  which 
method  would  be  better  in  what  situation.  But  this  would 
not  be  easy  to  accomplish.  As  Spronk  [S5]  has  noted, 
the  evaluation  of  the  appropriate  multiple  objective  model 
to  use  is  itself  a  multiple  objective  problem. 

The  purpose  of  this  section  is  not  to  give  those  general 
guidelines  in  the  selection  of  MODM  algorithms.  Instead,  it 
is  to  give  a  check  table  which  will  be  used  to  identify  the 
situations  under  which  GP  can  be  applied  correctly.  if 
GP  is  not  an  appropriate  method  to  be  used,  the  other 
methods  are  suggested.  In  the  following  subsections,  the 
different  types  of  MODM  methods  are  proposed,  the  criteria 
used  to  identify  the  situations  are  introduced,  and  finally 
a  choice  rule  for  GP  and  other  methods  is  presented. 
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4.2.1  The  Proposed  Alternate  Methods 

In  this  section,  some  MODM  methods  are  proposed  to 
replace  the  GP  method  in  those  situations  in  which  GP  is  not 
the  appropriate  method  to  be  used.  The  selection  of  these 
MODM  methods  is  based  on  their  representativeness  of 
different  types  of  methods.  These  methods  can  be  classified 
into  5  groups  according  to  the  stages  at  which  Information 
is  needed  and  the  types  of  the  needed  information.  No 
detailed  evaluation  has  been  done  in  the  selection  of  the 
proposed  methods.  The  purpose  of  selecting  those  different 
methods  is  only  to  show  that  GP  is  not  the  only  MODM  method 
and  should  not  be  used  in  all  situations.  There  are  many 
other  methods  which  are  more  suitable  than  GP  in  some 
situations.  The  principle  is  that  the  different  types  of 
problems  should  be  solved  by  different  types  of  methods. 

The  proposed  five  groups  of  methods  include  : 

(1)  Global  Criterion  Method  s 

This   method   does   not   need   any    interobjective  or 

subjective  preference   information   from   the   DM   once  the 

problem  constraints  and  objectives  have  been  defined.  Thus 

this  approach  requires  that  the  DM  be   able   to   accept  the 
solution  obtained  from  the  method. 

(2)  Utility  Function  Method  and  Bounded  Objective  Method  : 
This  group   of   methods   needs  the   cardinal   preference 

information  before  the  problem  can  be  solved.   The   decision 
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maker  must  give   some   judgement   about   specific   objective 
preference  levels  or  specific  trade-offs. 

(3)  Lexicographic  Method  and  Goal  Programming  : 

This  group  of  methods  also  need  the  preference 
Information  before  the  problem  can  be  solved.  But  the 
preference  Information  can  be  cardinal,  ordinal,  or  mixed. 
If  the  preference  information  is  mixed,  the  DM  must  rank  the 
objectives  in  the  order  of  their  importance. 

(4)  KSU-STEM  and  ISGP  II  :  (  Interactive  Methods  ) 

This  class  of  methods,  generally  referred  to  as 
"Interactive  Methods",  relies  on  the  progressive  definition 
of  the  DM's  preferences  along  with  the  exploration  of  the 
objective  space.  The  progressive  definition  takes  place 
through  a  DM-Machine  dialogue  at  each  iteration.  At  each 
such  dialogue,  the  DM  is  asked  about  some  trade-off  or 
preference  information  based  upon  the  current  solution  (  or 
the  set  of  current  solutions  )  in  order  to  determine  the 
final  compromise  solution. 

(5)  Multiple  Objective  LP  (MOLP)  and  Parametric  Method  : 
The   methods   of   this   group    determine    a    set    of 

nondominated  solutions  to  the  MODM  problem.  From  this 
subset  the  DM  choose  the  most  satisfactory  solution,  making 
implicit  trade-offs  between  objectives  based  upon  some 
previously  unindicated  or  nonquantif iable  criteria.   In   any 
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case,  the  trade-off  information  (  which  remains  implict  )  is 
received  from  the  DM  after  the  method  has  terminated  and  the 
subset  of  nondominated  solutions  has  been  generated.  This 
class  of  methods  does  not  require  any  assumption  or 
information  regarding  the  DM''s  utility  function. 

The  algorithms  of  the  methods  listed  above  can  be  found 
in  Hwang  and  Masud's  survey  of  MODM  methods  [H6]  [H7]  or 
other  reviews  of  multiobjective  techniques  [C6]  [G2]  [Rl] 
[Z3]. 
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4.2.2  The  Criteria  to  Identify  the  Situations 

The  suggestion  of  the  appropriate  methods  to  30lve  the 
MODM  problems  is  highly  dependent  upon  the  nature  of  the 
problems,  types  of  information,  stages  of  articulation  of 
preference  information,  DM's  familiarity  with  the  problems, 
and  DM's  preference  about  the  solution  process.  Several 
researchers  have  compared  various  techniques  along  various 
criteria,  including  Buchanan  and  Daellenbach  [B2],  Cohon  and 
Marks  [C5],  Rietveld  [Rl],  and  Wallenius  [W2].  For  example, 
the  criteria  used  in  Wallenius'  evaluation  of  three 
interactive  methods  [W2]  are  s 

(1)  DM's  confidence  in  the  best  compromise  solution 

(2)  Ease  of  use  of  the  method 

(3)  Ease  of  understanding  the  logic  of  the  method 

(4)  Usefulness  of  the  information  provided  to  aid  the  DM 

(5)  Rapidity  of  convergence,  measured  by  the  number  of 
cycles  and  the  total  time  for  solving  the  problem 

(6)  CPU-time  on  a  computer 

But  there  exists  no  clear  criteria  has  been  suggested  as 
a  unified  approach  in  selecting  the  MODM  methods.  These 
criteria  must  include  the  characteristics  of  diversified 
problems  and  methods,  different  decision  making 
environments,  and  DM's  ability  and  preferences.  The 
criteria  for  this  purpose  are  complex.  In  this  section,  the 
proposed  criteria  used  to  select  the  MODM  methods  are  mainly 
for  the  purpose  of  distinguishing  the   applicability   of   GP 
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from  that  of  other  methods.  These  criteria  are  listed 
below,  and  "value  function"  will  be  explained  later  in  the 
section. 

The  Proposed  Criteria  : 

(1)  What  is  the  type  of  value  function  of  the  objectives  ? 

(2)  Is  priority  information  of  objectives  available  ? 

(3)  In  which  3tage  is  the  preference  information   given  ? 

(4)  In  which  way  is  the  preference  modeled  ? 

(5)  What  is  the  number  of  solutions  presented  to  the  DM  ? 

(6)  Does  the  DM  prefer  interactive  approach  ? 


Value  Function  : 

The  "value  function"  has  the  same  meaning  as  the  utility 
function  except  that  value  function  is  used  only  when  the 
outcomes  are  certain.  So.  the  value  function  represents  the 
decision  maker's  values  about  objectives.  Hannan  [H2] 
suggested  that  the  development  of  value  function  should  be 
the  first  step  in  solving  the  MODM  problems.  The  following 
descriptions  will  follow  Hannan" s  explanation  of  value 
function. 

In  many  instances  the  Initial  objective  functions  are  in 
this  form  already,  but  frequently  they  are  not.  For 
example,  one  of  the  objective  functions  may  be  the 
nurse/patient  ratio  in  a  hospital,  and  the  decision  variables 
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may  be  the  number  of  each  type  of  nurse  to  employ.  The 
objective  function  i3  a  monotone  increasing  function  (the 
larger  the  decision  variables,  the  larger  the  objective 
function).  However,  the  corresponding  value  function  does 
not  always  increase  with  the  decision  variables.  Once  the 
nurse/patient  ratio  reaches  a  certain  point,  it  is  no  longer 
of  any  value  to  hire  more  nurses.  Thus,  a  value  function 
should  be  developed  to  generate  the  utility  of  hiring 
additional  nurses,  not  the  nurse/patient  ratio  as  a  function 
of  the  decision  variables.  "More"  is  not  always  "better" 
despite  the  fact  it  is  represented  that  way  in  the  objective 
function. 

The  following  discussion  describes  three  simple  types  of 
value  functions  as  shown  in  Figure  4.2.  In  each  of  the 
three  valuef unctions,  it  is  assumed  that  250  is  the  highest 
possible  value  for  the  corresponding  objective  function  and 
0  is  the  lowest  possible  value.  The  value  functions  attach 
a  utility  of  one  to  the  actual  goals  and  less  values  to  the 
points  further  away  from  the  goals. 

In  Value  Function  1,  the  decision  maker's  utility  is 
increasing  in  a  constant  manner  up  to  the  maximum  possible 
value  of  250.  If  GP  were  used  for  this  problem  by  choosing 
a  goal  less  than  250.  there  would  be  the  possibility  of 
obtaining  dominated  solutions.  This  type  of  problem  is  the 
most  controversial  application  of  GP. 

In  Value  Function  2.  the  decision  maker's  utility  is 
increasing  in  a  constant  manner  up  to  the   maximum   possible 
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Type  I 


Type  II 


Type  III 


Utility 


0        SO      100      150      200      2£0 


Objection 
Function  Vilue 


Utility 


i        Objection 

0         50     100      UD      200     250  Function  Velue 

Utility 


0         50      100      150 


1 Objection 

I  I 1       Functii 


200     250 


on  \'2lue 


Figure  '4.2  Example  value   function  for  an  objective 
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value  which  occurs  when  the  objective  function  is  100. 
After  this  point,  however,  all  numbers  have  the  same  maximum 
utility  of  1.  Here.  100  should  be  the  goal  and  only 
deviations  associated  with  values  less  than  100  should  be 
penalized.  For  this  example.  GP  is  an  appropriate 
methodology  to  use  In  conjunction  with  the  value  function. 

For  Value  Function  3.  The  value  150  should  be  the  goal, 
and  clearly  there  should  be  different  weights  assigned  to 
deviations  on  either  side  of  150.  Again,  GP  is  an 
appropriate  methodology. 

Obviously,  the  three  value  functions  given  above  are 
rather  simplistic  and  do  not  represent  all  possible 
real-world  preference  structures.  Nevertheless.  most 
complex  value  functions  can  be  characterized  as  one  of  three 
types  : 

(1)  monotone  increasing   everywhere   or   monotone   decreasing 
everywhere, 

(2)  monotone  increasing  followed   by  a  constant  function  or  a 
constant  function  followed  by  a   monotone   function,  and 

(3)  monotone  increasing  followed  by  monotone  decreasing   or 
vice  versa. 

The  process  of  developing  value  functions  can  make  the 
DM  more  familiar  with  the  problems.  For  the  objectives  that 
more  (less)  is  always  better  (  i.e.  the  value  functions  are 
monotonia  ),  GP  is  not  the   right   method.     If   the   value 
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functions  are  not  monotonic,  then  goals  other  than  the 
optimal  objective  function  values  are  preferable  and  GP  can 
be  used.  Adherence  to  these  rules  can  substantially 
simplify  the  determination  that  whether  the  GP  should  be 
used  or  not. 
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4.2.3   A  Choice  Rule  for  GP  and  Alternatives 

Baaed  on  the  proposed  methods  and  criteria  in  the 
previous  sections,  the  general  rules  of  selecting  the 
appropriate  methods  are  suggested  in  this  section.  Here, 
we  pay  more  attention  in  deciding  the  situations  in  which 
the  GP  is  applicable.  Table  4.3  shows  the  decision 
structure  for  selecting  the  MODM  methods .  Figure  4 . 3 
presents  the  same  decision  structure  by  a  hierarchical  flow 
chart.  The  proper  methods  are  found  by  answering  the 
subsequent  questions  : 

(1)  What  is  the  type  of  value  function  ? 

—  "type  1"  :  monotone  increasing  or  decreasing 

—  "type  2"  :  monotone  increasing  followed  by  a  constant 

function  or  vice  versa 
"type  3"  :  monotone  increasing  followed  by  monotone 
decreasing  or  vice  versa 

(2)  Are  the  goal  values  clear  ? 

"yes"  :  the  DM  can  clearly  specify  the  goal  values 
"no"   :  the  DM  can  not  specify  the  definite  goal 
values 

(3)  Does  the  DM  prefer  the  Interactive  solution   process   ? 

—  "yes"  :  interactive  methods  take  DM  longer  time  but 
have  higher  quality  of  final  solution 
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no"   :  non-interactive  methods  can  get  quick 
solution  but  it  may  not  satisfy  the  DM 

(4)  Is  priority  information  about  objectives  available  ? 

—  "yes"  :  the  DM  can  indicate  the  preference 

information  explicitly  or  implicitly 

--  "no"   !  the  DM  can  not  indicate  the  preference 
information 

(5)  Is  the  articulation  of  preference  information  a  posterior? 

"yes"  :  the   DM   can   not   supply   any   preference 

information  before  or  during  the  solution 
process  but  is  able  to  choose  his/her  best 
compromise  solution  from  a  set  of 
nondominated  solutions 
--  "no"  :  the  DM  can  indicate  his/her  preference 
before  or  during  the  solution  process 
explicitly  or  implicitly 

(6)  In  which  way  are  preferences  modeled  ? 

--  utility  function  (  relative  weights  among  objectives 

—  minimum  standards  (  for  each  objective  ) 

—  goal  values  (  for  each  objective  ) 

--  lexicographic  orders  (  of  objectives  ) 
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Table  4.3   Deci3ion  table  for  selecting  MODM  methods 


Situations 


IF    type  2  or  3  value  function  [a] 
AND   goal  values  are  clear  [b] 
AND  interactive  way  [c] 


Suggested  Methods 


ISGP  II    or 
KSU-STEM 


IF    [a]  AND  [b]  AND  not  Co] 
IF    [a]  and  not  [b] 


Goal  Programming 

ISGP  II    or 
KSU-STEM 


IF 
AND 

IF 

AND 

AND 


type  1  value  function  £d]  Global  Criterion 

no  priority  information  available  Method 

[d] 

priority  information  available  [e]  MOLP   or 

posteriori  articulation  of  preference  Parametric  Method 
information  [f] 


IF  £d]  AND  [e]  AND  not  ft] 

AND  utility  function  preference  [g] 

AND  [c] 

IF  [d]  AND  [e]  AND   not  [f] 

AND  [g]  AND  not  [c] 

IF  [d]  AND  [e]  AND  not  [f] 

AND  minimum  standards  [h] 

AND  [c] 


ISGP  II   or 
KSU-STEM 


Utility  Function 
Method 

ISGP  II   or 
KSU-STEM 


IF 
AND 

IF 

AND 

AND 

IF 
AND 

IF 

AND 

AND 

IF 
AND 


[d]  AND   [e]  AND   not  [f] 
[h]  AND   not  [c] 

[d]  AND  [e]  AND  not  [f] 

goal  levels  of  objectives  [i] 

[c] 

[d]  AND   [e]  AND   not  [f] 
[i]  AND   not  [c] 

[d]  AND  [e]  AND  not  [f] 
lexicographic  orders  [j] 
clear  goal  levels  [k] 

[d]  AND   [e]  AND   not  [f] 
[j]  AND   not  [k] 


Bounded  Objective 
Method 

ISGP  II   or 
KSU-STEM 


Goal  Programming 


ISGP  II   or 
KSU-STEM 


Goal  Programming 
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MODM  PROBLEM 


Value  Function 
Analysis 


Dl  : 
Type  of 
Function 


02  : 

Clear  goals  ? 


Q4  : 

Priority 
In-formation 
Available  ? 


Type  2  &  3 


Type 


Yes 


D5  : 

Posterior 
Articulation 
of  Preference 

"5 


No 


Interactive      ^^ 
way    ? 


No 


Global    Criterion 


MDLP 

Parametric  method 


Figure   4.3  Choice  rule  for  MODM  methods 
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Q6    :    In    which    way   are 
reference  modeled    ? 


Utility  -function 
(relative  weights) 


No 


D7    : 

Interactive 
way  ? 


ISGP   II 
KSU-STEM 


Utility 

function 

method 


2 

T 

Mini  mum 
standards 

No 

' 

QB  : 

Interacti v 
k.  way  ? 

Yes 

r 

1 

ISBP  II 
KSU-STEM 

' 

Bounded 

Objective 

method 

Figure   4.3  Choice  rule  for  MODM  methods    (cont'd.) 
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Q6  :  In  which  way  are 
reference  modeled  ? 


No 


Initial 
Goals 


09  : 
Interactive 

p.    way  ? 


ISBP   II 
KSU-STEM 


Goal 
Programming 


Lexicographic 
orders 


Yes 


BIO    : 

CI  ear 

p.   goals   ? 

No 


Lex i  cogr aphi  c 
method 


Goal 
Programming 


Figure   4.3  Choice   rule  for  MODM  nethods    (cont'd.) 
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The  above  analysis  shows  that  there  are  many  situations 
in  which  other  methods  are  more  suitable  than  GP.  GP 
has  very  limited  usage.  And  in  those  situations  in  which  GP 
is  applicable,  the  remedial  techniques  discussed  in  this 
chapter  are  still  suggested  to  be  used.  The  analysis  al30 
demonstrates  that  the  Interactive  methods,  i.e.  ISGP  II  or 
KSU-STEM.  can  be  used  in  many  different  situations.  We 
conclude  that  different  problems  should  be  solved  by 
different  methods.  The  analyst  should  be  very  careful  in 
deciding  the  appropriate  methods. 
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CHAPTER  5  Conclusions 

GP,  one  of  the  most  popular  MODM  methods,  has  been  a 
controversial  method  in  solving  the  MODM  problems.  In  the 
previous  chapters,  we  discussed  in  detail  the  advantages  and 
disadvantages  of  GP,  especially  the  disavdantages .  The 
purpose  is  to  clarify  the  inherient  shortcomings  of  GPand 
let  people  who  use  GP  understand  its  limitations  and  take 
some  necessary  steps  to  prevent  the  undesirable  results. 

The  GP  approach  relies  heavily  on  the  decision  maker's 
perceptions  of  the  range  of  choice  and  feasibility.  It 
requires  the  decision  maker  to  set  the  initial  goals  for 
each  objective.  But  in  most  practical  problems,  especially 
in  planning  problems,  the  value  a  goal  should  have  is  not 
apparent.  The  initial  goals  set  by  the  decision  maker 
usually  are  either  too  ambitious  or  too  conservative.  Since 
GP  is  based  on  the  satisficing  philosophy,  its  solutions  are 
very  sensitive  to  the  initial  goal  levels.  The 
inappropriate  initial  goals  often  result  in  a  dominated 
solution  or  nonpreferred  solution. 

The  preemptive  GP  disallows  tradeoffs  between  objectives 
of  different  priority  levels  and  is  proven  to  be 
incompatible  with  the  nonpreemptive  (  utility  preference  ) 
approach.  Some  nondominated  solutions  which  can  be  yielded 
by  the  nonpreemptive  approach  can  never  be  yielded  by  the 
preemptive  approach  with  any  ordering  of  the  goals.  The 
decision  maker  then  runs   the   risk   of   not   obtaining   the 
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preferred  solution.  Besides,  the  preemptive  algorithm 
greatly  reduces  the  range  of  acceptable  compromise  solutions 
by  sequentially  cutting  off  portions  of  the  feasible 
solution  space.  The  restrictive  nature  of  the  preemptive 
goals  may  totally  ignore  the  lower  priority  objectives. 
Thus,  preemptive  GP  should  be  used  only  when  the  priorities 
are  truly   preemptive. 

If  GP  is  the  desirable  technique  to  solve  a  MODM 
problem,  the  remedial  methods  proposed  in  this  study  can  be 
used.  The  PIS,  NIS,  MAR,  and  MAG  are  suggested  to  be  the 
reference  value  in  setting  the  initial  goals.  Two  different 
techniques  are  suggested  for  ensuring  the  nondominance  of  GP 
solutions.  The  nonpreemptive  approach  is  recommended  to 
avoid  the  shortcomings  of  preemptive  GP.  And  the 
eigenvector  method  is  suggested  to  simplify  the  weighting 
process:  If  possible,  the  interactive  GP  approach  would  be 
the  best  remedial  method  for  GP.  In  the  interactive  GP 
approach,  preemptive  weight  is  not  used,  the  setting  of  the 
initial  goal  is  also  not  so  objectionable,  and  the 
nondominated  GP  solution  is  guaranteed. 

Multiple  objective  decision  making  (  MODM  )  has  been  one 
of  the  fastest  growing  areas  of  OR/MS  during  the  last  15 
years.  Many  different  algorithms  have  been  developed  to 
solve  the  different  MODM  problems.  How  could  GP ,  with  the 
shortcomings  discussed  in  this  study,  still  be  so  popular  ? 
It  is  probably  because  the  shortcomings  of  GP  are  not 
well-known  and  GP  is  simple  and  easy  to   use.     A   lack   of 
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guidelines  in  helping  the  user  to  select  the  appropriate 
MODM  methods  may  be  another  reason.  The  choice  rule 
presented  in  this  study  is  a  rough  guideline  and  is  mainly 
used  to  identify  the  situations  in  which  GP  can  be  used. 
The  results  show  that  GP  should  be  used  only  in  limited 
situations.  Therefore,  the  user  of  GP  should  be  more 
careful.  From  a  managerial  point  of  view,  there  is  a  need  to 
have  an  evaluation  of  the  different  kind  of  methods  so  that 
the  user  can  determine  which  method  would  be  best  in  a  given 
situation.  Of  course,  this  would  not  be  easy  because  of  the 
diversified  MODM  problems  and  methods,  and  different 
decision  making  environments.  But,  this  will  have  to  be 
done  in  the  future  study  of  multiple  objective  decision 
making. 
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ABSTRACT 

Goal  Programming  (GP)  is  one  of  the  earliest  methods  in 
multiobjective  programming  and  has  been  applied  to  almost 
every  field  to  solve  the  MODM  problems.  As  an  extension  of 
Linear  Programming  (LP),  it  is  advantageous  for  its  simple 
formulation,  simple  algorithm,  and  easy  usage.  But  there 
are  some  Inherent  shortcomings  of  GP  which  are  not  apparent 
and  well-known  to  its  users. 

The  major  criticisms  of  GP  relate  to  the  use  of 
preemptive  priorities  and  the  a  priori  definition  of  goals. 
GP  requires  the  decision  maker  to  Indicate  the  initial  goals 
for  each  objective  before  solving  the  MODM  problem.  The 
criticism  concerns  the  difficulty  of  specifying  a  set  of 
goals  a  priori  without  more  knowledge  of  the  solution  space. 
With  regard  to  the  preemptive  priorities,  it  -is  contended 
that  the  use  of  preemptive  priorities  is  rigid  and 
inconsistent  with  utility  function  over  the  k  objectives, 
which  ultimately  leads  to  decisions  which  are  unacceptable 
to  the  decision  maker. 

Besides,  some  remedial  methods  are  suggested  to  enhance 
the  adaptability  of  GP.  And  efforts  are  also  made  to 
identify  the  conditions  in  which  GP  can  be  correctly 
applied.  If  GP  is  not  the  appropriate  method,  a  simple 
decision  rule  for  MODM  methods  is  proposed  to  recommend  the 
other  MODM  methods. 


